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Dark radiation and dark matter in supersymmetric axion models 
with high reheating temperature 
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Recent studies of the cosmic microwave background, large scale structure, and big bang nucle- 
osyntesis (BBN) show trends towards extra radiation. Within the framework of supersymmetric 
hadronic axion models, we explore two high-reheating-temperature scenarios that can explain con- 
sistently extra radiation and cold dark matter (CDM), with the latter residing either in gravitinos 
or in axions. In the gravitino CDM case, axions from decays of thermal saxions provide extra radi- 
ation already prior to BBN and decays of axinos with a cosmologically required TeV-scale mass can 
produce extra entropy. In the axion CDM case, cosmological constraints are respected with light 
eV-scale axinos and weak-scale gravitinos that decay into axions and axinos. These decays lead to 
late extra radiation which can coexist with the early contributions from saxion decays. Results of 
the Planck satellite will probe extra radiation at late times and thereby both scenarios. Further 
tests are the searches for axions at ADMX and for supersymmetric particles at the LHC. 

PACS numbers: 14.80. Va, 11.30.Pb, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 



Recent cosmological studies show trends towards a ra- 
diation content of the Universe at the onset of big bang 
nucleosynthesis (BBN) and much later that exceeds ex- 
pectations for standard three active neutrino species. 
The obtained limits on non-standard contributions AN e g 
to the effective number of light neutrino species N e g are 
still consistent with the standard value N e g ~ 3 at the 
1-2(7 level. However, BBN likelihood analyses based 
on recent studies of the mass fraction Y p of primordial 
helium [ll find posterior maxima of AN c g ~ 0.7- 
0.8 P,0,lj] and precision cosmology studies of the cosmic 
microwave background (CMB) and large scale structure 
(LSS) means of AN oS ~ 0.8-1.8 0-0 • While the BBN 
studies are limited by systematic errors (see e.g. @), the 
Planck satellite mission will provide soon results for N e g 
at the CMB decoupling epoch with an expected unprece- 
dented sensitivity of AiV e ff ~ 0.26 at the lcr level @, 0. 
In fact, the upcoming Planck results may turn the trend 
for extra radiation into evidence in the near future. 



Various explanations for AN c g 
plored in the literature invoking, e.g., 
nos I3I. [Toll . o ther light species [ll], [l2[ , 
tries 



1 have been ex- 
light sterile neutri- 
neutrino asy mmc- 
13|,[14|], or decays of heavy particles [41 1 la - 1281 ] . Here 



we study two classes of supersymmetric (SUSY) hadronic 
axion models which describe consistently extra radiation 
and cold dark matter (CDM) for a high reheating temper- 
ature after inflation of up to T R - 10 9 GeV or 10 11 GeV. 
In the considered (R-parity conserving) models, it may 
thereby be possible to generate the baryon asymmetry, 
e.g., via thermal leptogenesis with hierachical heavy Ma- 
jorana neutrinos [29(. Moreover, SUSY axion models are 
compelling since both the strong CP problem and the hi- 
erarchy problem are solved simultaneously. These models 
come with new fields including the axion a, the saxion a, 
the axino a, and the gravitino G, which can play impor- 
tant cosmological roles depending on their masses, the 



Peccei-Quinn (PQ) scale /pq, and the reheating tem- 
perature Tr. 

As the pseudo-Nambu-Goldstone boson associated 
with the U(1)pq symmetry broken spontaneously at 
/pq [13, IH|, the axion has interactions suppressed by /pq 
and a mass of m a ~ 6 meV(10 9 GeV//pq). With labo- 
ratory, astrophysical, and cosmological studies [H|, HH 
pointing to /pq > 6 x 10 s GeV, the axion is predicted 
to be an extremely weakly interacting particle (EWIP) 
with a tiny mass of m a < 10 meV. In SUSY set- 
tings, the saxion and the axino appear respectively as 
the scalar and the fermionic partner of the axion. They 
are EWIPs as well with masses m a and rria that de- 
pend on details of the model and of SUSY breaking. 
For example, one expects the saxion mass m a to be of 
the order of the gravitino mass mg in gravity- mediated 
SUSY breaking. As the gauge field associated with local 
SUSY transformations, the gravitino is another EWIP 
with interactions suppressed by the (reduced) Planck 
scale Mp = 2.4 x 10 18 GeV and a mass that depends 
on the SUSY breaking scale. While we do not assume 
a specific SUSY breaking model, m„ — mg is used in 
the main part of this work. Other than that, rag (to- 
gether with m a ) and ma are treated as free parameters 
set in a way to evade cosmological constraints. Model 
building aspects of the considered mass hierarchies will 
be considered elsewhere. 

In the first of the two classes that we consider, the 
gravitino is the lightest supersymmetric particle (LSP) 
that provides CDM. Here decays of thermal saxions into 
axions can provide AN e g ~ 0.8_prior to BBN 0, il 
ill, [H; see also [13, HIS, [II, 13, HI for extra radiation 
from late decays of non-thermal saxions. In the second 
class, a very light axino is the LSP, the gravitino the 
next-to-LSP (NLSP) and CDM resides in axions from 
the misalignment mechanism. Again, it is possible to 
have AiVeff ~ 0.8 from decays of thermal saxions into 
axions already prior to BBN. However, now there can be 
an additional contribution of AN e g ~ 1 but only well 
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after BBN from gravitino decays into the axion and the 
axino [l7],|2(j. F° r both classes, we show updated AN e g 
contours that will point to new limits on Tr once the 
Planck results on AN c g are available. Moreover, we de- 
vote particular attention to cosmological viability and to 
the interplay with present and potential future insights 
from SUSY searches at the LHC. 

Some points by which our present study goes beyond 
directly related existing studies 0, Ho| are the following. 
Decays are treated beyond the sudden-decay approxima- 
tion. In the G LSP case, the resulting AN e g contours 
are confronted explicitly with the Tr limit imposed by a 
gravitino density fig that cannot exceed the dark mat- 
ter density f2cDM- Here cosmological constraints require 
rria > 2 TeV such that axinos decay prior to the decou- 
pling of the lightest ordinary sparticle (LOSP), which de- 
notes the lightest sparticle within the minimal supersym- 
metric standard model (MSSM) . Axinos can then provide 
a sizable fraction of the total energy density of the Uni- 
verse when decaying and thereby produce entropy [36l — 
[38J. This is included in our calculations. Here we ap- 
ply an updated result for the axino abundance produced 
thermally in the early Universe, which we obtain by in- 
cluding quartic axino-squark-antisquark-gluino interac- 
tions [39] omitted in an earlier calculation [40j. In the 
a LSP case with the G NLSP, we present AN e g contours 
that account for both decays, G — > aa and a — > aa, ex- 
plicitly. Moreover, our treatment includes contributions 
of the gravitino-spin-3/2 components and of electroweak 
processes to the thermally produced gravitino yield. In 
both of the considered LSP cases, we account system- 
atically for the possibility that saxion decays into gluon 
pairs can have a sizable branching ratio and can thereby 
produce significant amounts of entropy. 

The remainder of this paper is organized as follows. 
In the next section we discuss the observational hints 
towards extra radiation beyond the SM and possible sce- 
narios depending on the upcoming Planck results. Sec- 
tion IIIII is devoted to general aspects of the considered 
SUSY hadronic axion models in high-X^ scenarios, which 
apply to the two explored LSP cases. This section con- 
tains our updated result for the primordial abundance of 
thermally produced axinos. The gravitino CDM and the 
axion CDM scenarios are presented in Sects. IIVI and [V] 
respectively. Here we consider the corresponding contri- 
butions to f2cDM, AN c ff, and entropy, provide resulting 
Tr limits, and address the testability of these scenarios. 
We summarize our conclusions in Sect. fVTl Appendix IA1 
provides details on our updated calculation of the ther- 
mally produced axino abundance, where hard thermal 
loop (HTL) resummation [4l|, |42| is used to treat screen- 
ing effects of the primordial plasma as in Ref. [4(|. In 
Appendix [B] approximate expressions for the numeri- 
cal results obtained in Sects. |IV] and [V] are given that 
allow for a qualitative understanding of those results. 
While m a = is assumed throughout the main part of 
this work, we briefly describe the changes that occur for 
m a 7^ itIq in Appendix [Cl 



II. EXTRA RADIATION 

One of our key motivations for the studies presented 
in this work is the trend towards extra radiation inferred 
from current cosmological investigations as summarized 
briefly in the Introduction. In this section we expand 
slightly on the description of the current situation and 
outline different possible perspectives depending on the 
outcome of the corresponding Planck results for AN c g. 

The standard model (SM) predictions of the total rel- 
ativistic energy density, 



1 + Mt 



(1) 



T at T ~ 1 MeV (be- 
e~ annihilation) and by 



are given by N c g = 3 and T„ = 
fore neutrino decoupling and e + 

Nes = 3.046 and T v = (4/ll) 1 / 3 T after neutrino decou- 
pling. Here p 1 is the photon energy density and Ti v ) the 
temperature of photons (neutrinos) . The effective num- 
ber of light neutrino species JV e ff increases slightly due to 
residual neutrino heating by e + e _ annihilation [43|. 

There are various ways to probe N e g and thereby non- 
standard contributions AN e g to which we refer as extra 
radiation. At the epoch of BBN, a speed-up of the Hub- 
ble expansion rate caused by AiV c g > leads to a more 
efficient 4 He output than in standard BBN. Observation- 
ally inferred limits on the primordial 4 He mass fraction 
Yp can thus be translated into AN e g limits. Much later, 
at the epoch of CMB decoupling, AN c g > affects the 
time of radiation-matter equality, leads to a less efficient 
early integrated Sachs- Wolfe effect, and reduces the scale 
of the sound horizon. This affects the CMB power spec- 
trum by increasing the height of the first peak and by 
shifting the peak positions towards higher multipole mo- 
menta. Moreover, free-streaming of the relativistic pop- 
ulations associated with AN e g > supresses power on 
small scales and thereby affects the matter power spec- 
trum inferred from studies of the LSS. Based on those ob- 
servables, numerous studies of BBN, CMB, and LSS have 
explored limits and favored values for AJV e ff 0-0, HiHIH 
with the outcome outlined in the Introduction. 

To motivate the AN e g values considered in our study, 
we quote representative current constraints on AN e g im- 
posed by BBN and precision cosmology in Table |U The 
first two lines have been obtained in a BBN-likelihood 
analysis 4] based on the recent Yp studies of Izotov and 
Thuan [Jand of Aver et al. @ ■ Those studies report pri- 
mordial 4 He abundances of Yp IT = 0.2565±0.001(stat.)± 
0.005(syst.) and Y p Av = 0.2561 ± 0.0108, respectively 
with errors refering to 68% intervals. Moreover, a pri- 
mordial D abundance of log[D/H] p = -4.56 ± 0.04 |47j 
and a free- neutron lifetime of r n = 880.1 ± 1.1s [33[ have 
been used in the determination of the listed posterior 
maxima (p.m.) and the 3ct upper limits. The third line 
gives the mean and the 95% confidence level (CL) upper 
limit on AN e g obtained in the precision cosmology study 
of Ref. [6] based on CMB data, the Sloan Digital Sky Sur- 
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TABLE I. Constraints on AN c g from BBN and precision cos- 
mology. The first two lines give the posterior maximum (p.m.) 
and the minimal 99.7% credible interval imposed by BBN as 
obtained in Ref. [J| using the indicated data sets and the prior 
AiVeff > 0. The third line lists the mean and the 95% CL up- 
per limit on AA^ff from the precision cosmology study [y| 
based on CMB data, the Sloan Digital Sky Survey (SDSS) 
data-release 7 halo power spectrum (HPS), and data from 
the Hubble Space Telescope (HST). 



Data 


p.m. /mean 


upper limit 


F P IT [J] + [D/H] p [47] 


0.76 


< 1.97 (3a) 


Y p Av \2] + [D/H] p [47] 


0.77 


< 3.53 (3a) 


CMB + HPS + HST [6J 


1.73 


< 3.59 (2a) 



vey (SDSS) data-release 7 halo power spectrum (HPS), 
and data from the Hubble Space Telescope (HST). Com- 
paring this mean with the posterior maxima quoted in 
the first two lines, one sees that the favored AN c g value 
at the CMB decoupling epoch exceeds the one at the 
BBN epoch by about one. Moreover, the conservative 3cr 
BBN-imposed upper limit is very similar to the 2a up- 
per limit from precision cosmology. Compatibility with 
AN cS = is found at the l-2a level in both the BBN 
and the precision cosmology study. While a more de- 
cisive compatibility test seems to be difficult for BBN 
investigations due to significant systematic uncertainties 
(see e.g. [2j), the results of the Planck satellite mission 
are expected to improve the AN e g accuracy of precision 
cosmology investigations substantially. As mentioned in 
the Introduction, a la sensitivity of AN e g ~ 0.26 has 
been found in error forecasts for future precision cosmol- 
ogy studies [1, . This may clarify the extra-radiation 
situation at the CMB-decoupling epoch. 

If it remains unclear whether an excess of radiation 
was already present at the BBN epoch, there will be the 
following possibilities depending on the outcome of the 
Planck satellite mission: 

• Planck finds no statistically significant excess, i.e., 
AN c g < 0.78 at the 3a level. Nevertheless, there 
can still be a small amount of extra radiation, e.g., 
from saxion decays and/or gravitino decays. With 
the current BBN limits, the following scenarios are 
possible: (i) this small amount was already present 
at the onset of BBN with no additional contribu- 
tion after BBN, (ii) this small amount was gener- 
ated only well after BBN, or (iii) part of this small 
amount was generated already prior to BBN and 
the remaining part well after BBN. 

• The ambiguous situation with only a hint towards 
extra radiation will remain if the new 3a limit is 
AN c g < 1.56. Again compositions (i)-(iii) but now 
of the possibly larger AN e g are possible. 

• Planck finds a statistically significant excess of dark 
radiation at the 3a level, i.e., a mean larger than 



AiVeff ~ 0.78 and a 3a limit of AN e g > 1.56. Ac- 
cordingly, a significant amount of extra dark ra- 
diation can be provided by saxion decays and/or 
gravitino decays. Once more the three composi- 
tions (i)-(iii) described above will be possible. 

We will see below that composition (i) is the only 
one that can be realized in the considered gravitino LSP 
case, whereas the alternative axino LSP case allows for 
all three compositions. Contours of AN c g — 0.26, 0.52, 
0.78, 1.73, and 3.59 will be explored in the respective 
parameter regions corresponding to the forecasted la, 
2a and 3a Planck accuracies and to the mean and 2a 
limit of the precision cosmology study quoted in Table HI 
Accidentally, the value 0.78 indicating 3a Planck accu- 
racy nearly coincides with the posterior maxima from 
the BBN analysis quoted in Table |TJ Thus, we will use 
the respective contours to infer also parameter regions in 
which one finds the AN e g value favored by BBN studies. 
Moreover, we will consider contours of AN c g = 1.73 also 
as an exemplary hypothetical 3a upper limit to illustrate 
the possible case in which Planck results point to extra 
radiation in excess of the statistically significant 3a level. 



III. HIGH-REHEATING-TEMPERATURE 
SCENARIOS 

Throughout this work it is assumed that inflation has 
governed the earliest moments of the Universe, as sug- 
gested by its flatness, isotropy, and homogeneity. Ac- 
cordingly, any initial EWIP population was diluted away 
by the exponential expansion during the slow-roll phase 
of the inflaton field. A radiation-dominated epoch with 
an initial temperature of Tr emerged from the subse- 
quent reheating phase in which inflaton decays repop- 
ulate the UniverseQ While inflation models may point 
to Tr well above 10 10 GeV, we limit our studies to the 
case Tr < /pq in which no PQ symmetry restoration 
takes place after inflation. Focussing on high-reheating 
temperature scenarios with Tr > 10 7 GeV, axions, sax- 
ions, axinos, and gravitinos can be produced efficiently 
in thermal scattering of MSSM fields in the hot plasma. 
Depending on the PQ scale /pq and on Tr, even scenar- 
ios in which the fields of the axion supermultiplet were 
in thermal equilibrium are conceivable. 

For the axion and the saxion, our estimate for the de- 
coupling temperature reads [|| 



T£ a « 1.4 x 



10 9 GeV 



/PQ 



10 11 GeV 



(2) 



Following the approach of Ref. [4j and using our results 
for thermal axino production presented below and in Ap- 



1 Inflaton decays into EWIPs may have been efficient. However, 
we do include such contributions since there are inflation models 
in which this production mechanism can be negligible 48, 49]. 
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pendix^l we estimate the axino decoupling temperature 
as 



rpa 

1 U 



5.2 x 10 8 GeV 



fpQ 



10 n GcV 



(3) 



In cosmological scenarios with Tr > T-g (or even 
Tr > Tp' CT ), axinos (together with axions/saxions) were 
in thermal equilibrium before decoupling as a relativistic 
species provided ma <C (and m a C 1^). Then the 
yield of those thermal relic axions/saxions and axinos af- 
ter decoupling is given respectively by 



TABLE II. Assignments of the index i, the gauge coupling 
gi, and the gaugino mass parameter Mi, to the gauge groups 
U(1)y, SU(2)l, and SU(3) C , and the constants ki, yi, and tJi. 



gauge group 


i 


9i 


Mi 


ki 


(y*/i(n 14 ) 




U(1) Y 


1 


9 


Mi 


1.266 


0.653 


0.018 


SU(2) L 


2 


9 


M 2 


1.312 


1.604 


0.044 


SU(3) C 


3 


9s 


M 3 


1.271 


4.276 


0.117 



,eq 



Y eq = 



1.2 x 10" 



and 



Y, 



cq _ 



1.8 x 10 



-3 



(4) 



(5) 



Here denotes the corresponding (equilibrium) num- 
ber density of species j and s the entropy density. For the 
latter, we use s(T) = 27r 2 <7*sT 3 /45 with an effective num- 
ber of relativistic degrees of freedom of g*s(Ib) — 232.5 
that accounts for the MSSM and the axion multiplct 
fields, which can all be considered as relativistic at Zb 
for m a ,a < T D . 

In scenarios with Tr, < T^ a ' a , the axion multiplet 
fields can still be thermally produced (TP) via scatter- 
ing of colored (s)particles in the primordial plasma. The 
resulting yields are given by [J] 



F a T P = 1.33x 10- a ff «ln 



L01 \ / 10 u GeV \Y Tr 
~J\ /pq J V10 8 GeV 
_(6) 

and, as derived by updating the result of Ref. [40j in 
Appendix El by 



Y? v = 1.98xl0^1n 



1.27 V 10 11 GeV Y 



9s 



fpQ 



10 s GeV 



(7) 

Here the strong gauge coupling is understood to be eval- 
uated at T R , i.e., g s = g s (T R ) = v /47ra s (T R ), which we 
calculate according to its 1-loop renormalization group 
running within the MSSM from a s (mz) = 0.1176 at the 
Z-boson mass mz = 91.1876 GeV. 

Note that our focus is on hadronic or KSVZ axion mod- 
els (ESIEl in a SUSY setting, _with N Q = 1 heavy KSVZ 
(s)quark multiplcts Ql and Qr. After integrating out 
the KSVZ fields, we obtain the effective Lagrangian 0] 



where b is a color index, the corresponding color- 
gauge covariant derivative, G h the gluon-field-strength 

tensor, G b — e l _ Ll/p<J G b pa /2 its dual, g b the gluino field, 
and D b = — g s J^- with a sum over all squark 

fields q and the SU(3) C generators Th in their fun- 
damental representation; the subscript M indicates 4- 
component Majorana spinorsQ In the considered frame- 
work, the Lagrangian ((8|) describes the relevant sax- 
ion/axion/axino interactions even in a conceivable very 
hot early stage of the primordial plasma with tempera- 
tures T not too far below /pq0 Based on © the pre- 
sented results ©, ©, ©, and © are obtained. In 
particular, as outlined in more detail in Appendix [XI 
our result for the thermally produced axino yield (© ac- 
counts for the second term in the third line of © that de- 
scribes the quartic axino-squark-antisquark-gluino inter- 
action 39] , whereas the corresponding result of Ref. [4(| 
was based on only the first term in that line. 

Gravitinos with mass values of mg > 1 GeV, which 
are the ones considered in this work, have never been in 
thermal equilibrium with the primordial plasma. Never- 
theless, they can be produced efficiently in thermal scat- 
tering of MSSM fields in the hot plasma. Derived in a 
gauge- invariant treatment, the resulting thermally pro- 
duced gravitino yield reads [H, H?! 



(9) 

with yi, the gauge couplings gi, the gaugino mass param- 
eters Mi, and ki as given in Table HT1 Here Mi and gi are 
understood to be evaluated at Tr. 

In the following we consider universal gaugino masses, 
7774/2 — Mi(mGUT), at the grand unification scale 
777GUT ~2x 10 16 GeV. We do not specify a SUSY model. 




*-PQ 



a {G b ^G b - 2D b D» - 2tfVy*£> M &) 



^ + 2g° M ^7°D fl g b M 



a\G b ^G b -u •>/,':.-/'- 



ia M - 



2a mD g M 



(8) 



Slightly different expressions for £J?q can be found in [39l . [52 



We use the space-time metric </ M „ = 



■■ diagf+1, -1, -1, -1) 



and other conventions and notations of Ref. 1531 and, except for 
a different sign of the Levi- Civita tensor e 0123 = +1, of Ref. @. 
3 We do not consider scenarios with a radiation-dominated epoch 
with T above the masses of the heavy KSVZ (s)quarks rrig q 
such as those considered in Ref. |55||. 
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Nevertheless, we use certain pairs of and the weak- 
scale gluino mass rrig keeping in mind that these values 
are related via renormalization group evolution. In par- 
ticular, we will associate nig ~ 1, 1.25, and 1.5 TeV with 
m i/2 — 400, 500, and 600 GeV, respectively. Computing 
the renormalization g rou p evolution with the spectrum 
generator SPHENO [53, [59|, these relations are obtained 
within the Constrained MSSM (CMSSM) with a univer- 
sal scalar mass parameter of mo = 1.7 TeV, the trilin- 
ear coupling A Q — 0, a positive higgsino mass parame- 
ter, /i > 0, and a mixing angle in the Higgs sector of 
tan/? =10. The above combinations are still allowed by 
current SUSY searches at the LHC but are well within 
reach of the ongoing experiments; see e.g. Re f. [ 60l| . 

Note that the field-theoretical methods [HI, |42| applied 
in the derivations of ([5]), ©, and ^ require weak cou- 
plings Qi < 1 and thus T > 10 6 GeVQ Moreover, in 
those derivations, a hot thermal plasma consisting of the 
particle content of the MSSM is considered in the high- 
temperature limit. In fact, it is assumed that radiation 
governs the energy density of the Universe as long as 
thermal production of the respective EWIP is efficient, 
i.e., for T down to at least T ~ 0.01 Tr. This is assumed 
in this work also. However, we will encounter situations 
with significant entropy production at smaller temper- 
atures generated by decays of by then non-relativistic 
saxions and/or axinos from thermal processes. Then this 
can dilute the yield of a stable or long-lived EWIP from 
thermal processes in the earliest epoch correspondingly 
with dilution factors of A > 1: 



oq/TP 



EWIP 



1 

A 



cq/TP 



EWIP 



(10) 



Abundances of decoupled species that emerge from de- 
cays of thermally produced EWIPs prior to the entropy 
producing event are equally affected. 

In high-reheating temperature scenarios, the LOSP 
usually freezes-out as a WIMP at a decoupling tempera- 
ture of Tp OSP ~ m L osp/25 with an abundance Flosp 
that can be determined by solving the corresponding 
Boltzmann equations. In the case of entropy production 
after LOSP decoupling, this abundance will be diluted 



Vlosp — > ^-Vlosp 



(11) 



as well. However, in situations in which the entropy pro- 
ducing event ends well before LOSP decoupling, Ylosp 
is not affected. Here we assume in both cases that LOSP 
decoupling takes place in a radiation-dominated epoch. 
This is justified in the settings considered below where 
the contribution of long-lived non-relativistic species to 
the total energy density (that enters the Friedmann equa- 
tion) is negligible during LOSP freeze-out. 

In high-reheating temperature scenarios, thermal lep- 
togenesis with hierarchical heavy Majorana neutrinos 



can explain the baryon asymmetry of the Universe [29|. 
Without late-time entropy production, Mri ~ Tr of 
at least about 10 9 GeV is then required to generate the 
observed baryon asymmetry rj, where Mr\ denotes the 
mass of the lightest among the heavy right-handed Ma- 
jorana neutrinos. With late-time entropy production, 
a baryon asymmetry generated prior to the entropy- 
producing events must have been larger by the associ- 
ated dilution factor A. In the framework of thermal lep- 
togenesis, this can be realized for up to A ~ 10 4 with 
Mr-i ~ T R . ~ 10 13 GeV, as can be seen in Fig. 7(a) of 



; see also 



In 



Ref. [62j and in Fig. 2 of Ref. 
fact, with a dilution factor of A, the required minimum 
temperature for successful leptogenesis has to be larger 
by that factor: 



Tr > 10 9 GeV 



1 

A 



Tr > 10 9 GeV. (12) 



Together with (flO|) and (fTTj) . this motivates us to care- 
fully calculate A and to monitor the results for the two 
scenarios discussed in the following. 



IV. GRAVITINO CDM CASE 

In this section we look at the R-parity conserving 
SUSY scenario in which a gravitino with mass mg > 
1 GeV is the stable LSP whose thermally produced den- 
sity parameter 

nfh 2 = m d Y? p (T )s(T )h 2 /p c (13) 

provides the dominant part of the CDM density [33j 

ficDM/i 2 = 0.111 ±0.018, (14) 

with errors at the 3cr level, the present photon temper- 
ature To = 0.235 meV, the Hubble constant h in units 
of 100 kmMpc-V 1 , and p c /[s(T )h 2 } = 3.6 eV. Ac- 
cordingly, all heavier sparticles including the LOSP and 
the axino are unstable. In turn, each LOSP and each 
axino present in the Universe after LOSP decoupling 
will decay directly or via a cascade into one gravitino. 
Depending on Ilosp? the contribution to f2g from de- 
cays of thermal relic LOSPs can be small as will be dis- 
cussed below in more detail. This is different for long- 
lived axinos that decay at temperatures below a fidu- 
cial T\ ow <C Tp OSP . For the considered settings with 
fff p ~ S! CDM and /pq < 10 12 GeV, their contribution 



(Ti ow )s(T )h 2 /p c 



(15) 



For alternative approaches, sec 39, 52, 61 



exceeds (|14|) by many orders of magnitude. This can be 
immediately seen when comparing <j5j) and © with 
To avoid this excess, we focus in this section on G LSP 
scenarios in which axinos decay dominantly into gluons 
and gluinos well before LOSP decoupling with a rate that 
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can be derived from the effective Lagrangian 



T~ ~ v. 



■U'l 



167T3/| Q 



1-4 . a*) 



In the G LSP scenarios considered in this section, ax- 
ions from decays of thermal saxions prior to BBN are 
the only contribution to A.N e s, as already mentioned in 
Sects. UlandHIl The Lagrangian that allows for the rele- 
vant a — > aa decay reads |15| 



/•kin 

^PQ — 



\/2x \ 
+ a 

VPQ J 

2 2 



and the associated decay rate 



(17) 



32ntf 



(18) 



PQ 



where x — £V qfvf /vpq depends on the axion model with 
qi denoting the charges and Vi the vevs of the fundamen- 
tal PQ fields [15(. For example, x = 1 in a KSVZ axion 
model with just one PQ scalar (with q = 1 and v — i>pq) 
and x <C 1 in such a model with two PQ scalars with 
qi = —q2 = 1 and similar vevs, v\ ~ v-i ~ wpq/v^- 
The two scales zjpq = \/^vfq^ and /pq are related via 
f P Q = V2vpq [4] . 

For m a > 1 GeV, the saxion decay into two gluons, 
a — > gg, can become a competing decay mode towards 
small values of x. The associated rate reads 



°^gg 



(19) 



and is derived from ((5J) . The saxion decay into photons is 
subdominant whenever the a — > gg decay is kinematically 
viable, i.e., for m a above the threshold to form hadrons. 
Saxion decays into gluinos or axinos are kinematically not 
possible in the G LSP case with m a = tuq. Accordingly, 
the lifetime of the saxion and the branching ratio of its 
decays into axions and into gluons are well described by 



1 



1 



327T/2 



PQ 



BR(er — S- aa) 



X 



°^gg 

2 



ml[x 2 +2(a s /7r) 2 



+ 2(a s /7r) 2 ' 
BR(er -> gg) ~ 1 - BR(a -> aa) = 



1 



2(a s /^) 2 



(20) 
(21) 

n /(22) 



respectively, with a s = a s (r7i cr ). For example, for x > 0.2 
and vn a > 10 GeV, one finds BR(ct — > aa) > 0.9 so that 
To- is governed by the decay into axions. Towards smaller 
x and/or m CT , the saxion decay into gluon pairs becomes 
important with effects discussed below. 



When decaying, both the axino and the saxion are non- 
rclativistic. Accordingly, we encounter two types of de- 
cays of non-relativistic particles: (i) decays into axions 
which are by then decoupled from the thermal plasma 
and thereby an inert relativistic species and (ii) decays 
into relativistic species that are rapidly thermalized and 
thereby associated with entropy production. We can in- 
deed face simultaneously situations studied previously for 
the generic cases of out-of-equilibrium decays of non- 
relativistic particles into inert radiation [64j and into 
thermalizing radiation that produce entropy [65] . 

Let us now calculate the contribution to AN e g of the 
energy density of relativistic axions p a from thermal pro- 
cesses in the earliest moments of the Universe and - more 
importantly - from late decays of thermal saxions, 



AAf eff (T) 



120 



lPa{T), 



(23) 



and the relic density of thermally produced gravitinos 
ilj^h 2 . By taking into account the possibility of entropy 
production in both axino and saxion decays, we general- 
ize and refine our related previous study Q. Moreover, 
our numerical results are now obtained beyond the sud- 
den decay approximation. Nevertheless, we will return 
to that approximation to derive expressions that allow 
for a qualitative understanding of the behavior of our 
numerical solutions in Appendix [Bj 

In the epoch when thermal processes involving EWIPs 
are no longer efficient and when axinos and saxions from 
such processes are non-relativistic, the time evolution of 
the energy densities of axinos, saxions, and relativistic 
axions is described by the following Boltzmann equations 

f>a+3H P a = -TaPa, (24) 
p a + "ill Per = -Y a p a , (25) 
p a + AHp a = BR((T -> aa)T a Pa, (26) 

with the Hubble expansion rate H = R/R and the dot 
indicating derivation with respect to cosmic time t. The 
time evolution of entropy S is given by 

X {TaPa + [1 - BR(<7 -+ aa)]T aP a} , (27) 

and the one of the cosmic scale factor R by the Friedmann 
equation for a flat Universe 



H 2 



3to p 



(pa + Pa + Pa+ Prad) 



(28) 



with the Planck mass rap = 1.22 x 10 19 GeV and the en- 
ergy density of the thermal MSSM radiation background 



3 3* 



30 4g tS \2TT z g*s 



45 



1/3 



SV3 

R 4 



(29) 



where g* is the effective number of relativistic degrees of 
freedom within the MSSM only, i.e., without the axion 
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multiplet and the gravitino. In this section, = g*s 
holds for the interval over which we integrate the Boltz- 
mann equations. (This will be different in Sect EI ) 

Equations (fM)) ~ (f2"8|) form a closed set of differential 
equations that we solve numerically. We begin our com- 
putation at t{ — 1.6 x 10~ 13 s corresponding to T = 
1 TeV with R{U) = 1 GeV -1 and end at tf = 0.7 s cor- 
responding to Tf ~ 1 MeV. For the initial values of the 
energy densities, we use 

PaiU) = m a rf q/TP S (T t ), (30) 
p a (U) = m^/^siTi), (31) 
Pa(U) = (pZ)Ya q/TP s(Ti), (32) 

where the average thermal axion momentum is (p„ h j) = 
2.701 T a>i and T a>i = [g., s {T. l )/228.7b] 1 / 3 T i . Note that 
saxions can be treated as a non-relativistic species 
throughout the time interval [ti, tf] although a saxion, 
e.g., with Too- = 100 GeV will be relativistic at an initial 
temperature of T — 1 TeV. At times at which saxions 
are relativistic, their contribution p a on the right-hand 
side of the Friedmann equation (l28|) is negligible. When- 
ever their contribution becomes sizable, they are non- 
relativistic, which justifies the simplified treatment. 

With the initial entropy S(ti) = s(Ti)R(ti) 3 and after 
numerical integration, we obtain the dilution factor 

As described already in the previous section, this factor 
quantifies the dilution due to entropy release which af- 
fects the yield of species not in thermal equilibrium such 
as Y~ p and thereby (|13p . In fact, since p~ p can be safely 
neglected in (|28|) at the considered times and since the 
gravitino is stable in the case considered here, it is not 
necessary to include the Boltzmann equation for the grav- 
itino in the described calculation. This will be different 
in the axino LSP case considered in the next section. 

Results of our numerical integration are illustrated in 
Fig. CD for m~ a = 6 TeV, m~ g = 1 TeV, and / PQ = 
10 11 GeV. For this setting, T aftGr ~ 10 GeV is the tem- 
perature at the end of the axino-decay epoch, at which 
Ta — 2>H is satisfied. Therefore, a realistic LOSP with 
tolosp ^ 250 GeV is compatible with the requirement 
Tkftor > Td OSP (~ to L osp/25) that is crucial as dis- 
cussed at the beginning of this section. As mentioned 
in Sect. IIII1 the considered gluino mass is still compati- 
ble with limits from SUSY searches at the LHC [60j . 

Figure[T](a) shows the time evolution of the energy per 
comoving volume, R 3 p, of axinos (dash-dotted), saxions 
(dashed) , axions (dotted) and other radiation (solid) and 
of entropy S (dash-double-dotted) for m a = 100 GeV 
and Tr = 10 9 GeV, where black and gray lines are ob- 
tained with x = 1 and 0.02, respectively. In both cases 
there is extra radiation in the form of axions. Consid- 
ering the dotted lines prior to saxion decay, one can see 



clearly that only a very minor contribution resides in ax- 
ions from thermal processes. Thus, the extra radiation 
AA^ff resides basically fully in axions from saxion de- 
cays, as can be seen by the rise of the dotted lines that 
results from those decays. For x — 1, one can see that 
entropy with a dilution factor of A ~ 2 is produced in 
axino decays only. This is different for x = 0.02 where 
additional significant late contributions to S and thereby 
to A emerge from saxion decays into gluons. Consider- 
ing i? 3 p ra d, one sees that it decreases slower during the 
entropy producing event (s), whereas other dips of that 
solid line result from changes in the effective number of 
relativistic degrees of freedom. In contrast to the ax- 
ino, which is required to decay prior to LOSP freezeout, 
entropy released in late saxion decays dilutes Ilosp m 
addition to, e.g., Y~ F or the baryon asymmetry. 

Or 

In general, towards small x, both the saxion lifetime T a 
and BR(ct — > gg) increase which leads to larger values of 
A. This effect becomes even more pronounced towards 
smaller m CT as long as the decay a — > gg is not kincmat- 
ically suppressed. Figure [TJb) illustrates this behavior, 
which also becomes manifest in the approximations (|B6[) 
and (|B7[) obtained in Appendix [B] This panel of Fig. []] 
shows the dilution factor A as a function of the reheat- 
ing temperature Tr. Here black and gray lines refer to 
m a = 20 and 100 GeV, respectively, and are presented for 
x = 1 (solid), 0.2 (dashed), 0.1 (dotted), and 0.02 (dash- 
dotted). The Tr dependence of A results from the one of 
V a cq/TP and of V CT oq/TP ; cf. {B5|) - (|B7j. in Appendix!] The 
kink in the A contour that is visible already for x = 1 
indicates the Tr value that coincides with the decoupling 
temperature of axinos T^ given in ([3]). In cosmological 
scenarios with Tr > Tg, Y~ q applies which is indepen- 
dent of Tr. The other kinks at larger Tr that appear for 
a;<l indicate the corresponding Tr value above which 
Tr > Tp, where the latter is given in @. With a more 
careful treatment that includes axino/saxion disappear- 
ance reactions when calculating the thermally produced 
yields for Tr near the respective decoupling tempera- 
tures, these kinks will disappear. Expecting smoother 
curves that are close the shown ones, we leave such a 
treatment for future work. 

Let us now explore systematically the amount of extra 
radiation released by saxion decays and regions in which 
the constraint f2~ P < ficDM is respected. Results for 
x = 1 are presented in Fig. [2] and for x = 0.1 and 0.2 
in Fig. [3] In both figures, we consider m a = tUq and 
m±/2 = Mi (tigut)' As already discussed in Sect. [TTJ 
there are hints towards the existence of extra radiation. 
These hints could be an indication for the existence of 
axions from saxion decay. We investigate this possibil- 
ity for f PQ = 10 10 , 5 x 10 10 , and 10 11 GeV. For each of 
these values, ma and mg are chosen such that the ax- 
ino decay can take place before the freezeout of a not too 
massive LOSP. We report the considered combinations in 
Table IIIII together with T a f tG r at which Ta = 3H and the 
mass of a LOSP to™osp f° r which its decoupling temper- 
ature satisfies T^ OSP ~ to LO sp/25 = T after . This table 
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FIG. 1. (a) Time evolution of the energy per comoving volume, R 3 p, of axinos (dash-dotted), saxions (dashed), axions (dotted) 
and other radiation (solid) and of entropy S (dash-double-dotted). Here m a = 100 GeV, ma = 6 TeV, nig = 1 TeV, 
Tr = 10 9 GeV, and /pq = 10 11 GeV. The initial value of the scale factor is set to Ri = 1 GeV -1 at the temperature Tj = 1 TeV 
corresponding to a time of ti = 1.6 x 10~~ 13 s. Black (gray) lines refer to the case with x = 1 (0.02). (b) The dilution factor A 
as a function of the reheating temperature Tr for x = 1, 0.2, 0.1, and 0.02 shown by the solid, dashed, dotted, and dash-dotted 
lines, respectively. Black (gray) lines are obtained with m CT = 20 (100) GeV, whereas all other parameter are as in panel (a). 



TABLE III. The temperature T aftor at which F a ~ 3H for 
different combinations of the PQ scale /pq, the axino mass 
ma, and the gluino mass nig together with the LOSP mass 
m5§|p for which T£ OSP ~ m LOS p/25 ~ T after . 



/pQ 


nia 


m g 


2~aftcr 


m max 
m LOSP 


[GeV] 


[TeV] 


[TeV] 


[GeV] 


[GeV] 


10 10 


2 


1 (1.25) 


13 (9) 


325 (225) 


5 x 10 10 


3 


1 (1.25) 


6 (5) 


150 (135) 


10 11 


6 


1 (1.25) 


10 (9) 


250 (235) 



shows explicitly that the viability of these gravitino LSP 
scenarios requires the axino to be quite heavy and the 
LOSP to be relatively light. 

Figures [2Ia)-(c) show the amount of extra radiation 
AiV e ff provided by axions from decays of thermal sax- 
ions for x = 1 together with the upper limit on Tr 
imposed by fl]~ p h 2 < 0.129 at the 3tr level. The solid 

Cr 

black (gray) contours indicate J7~ p /i 2 = 0.129 for mi/2 = 
400 (500) GeV. The high Tr regions above these contours 
are disfavored by overly efficient gravitino production. 
The dashed, dotted, and dash-dotted contours indicate 
AiV cfi = 0.78, 0.52, and 0.26, respectively. Here black 
(gray) contours are obtained with m,g = 1 (1.25) TeV, 



which is compatible with mi/ 2 = 400 (500) GeV used to 
evaluate ^i~ p . The Tr dependence of the AiV e ff contours 
disappears for cosmological scenarios with Tr > T^. The 
difference between the black and gray AN e g contours for 
a fixed AN c g results from the dependence of the dilution 
factor A on uig. The corresponding dilution factors A 
can be read from Fig. Hfd). The enhanced kinematical 
suppression of axino decays for a heavier gluino leads 
to a longer axino lifetime and thereby to a larger A, 
which also can be seen in (|B5|) of Appendix [B] This di- 
lutes Ya q ^ TF or p a more strongly and thus reduces AN e g 
correspondingly at a given combination of m a and Tr ; 
cf. dB8j) and (jBTOf in Appendix|U Moreover, in Figs.^b) 
and (c), one can see kinks in the AiV e ff contours at Tr, 
values below T^. Those kinks appear at the same Tr 
values as the kinks in Fig. [U[d) and indicate the point 
above which Tr > Tg. The dilution factors obtained for 
nig = 1 and 1.25 TeV are also included in the calculation 
of the f2~ p contours for m 1 / 2 — 400 and 500 GeV, respec- 

tively. Indeed, the kinks in the i7~ p contours visible in 
the panels (b) and (c) result from the A behavior shown 
in panel (d). Despite the larger A for larger mg, the Tr 
limit imposed by f2~ p < Hcdm is still more restrictive 
for larger due to the Mj dependence of ©. 

As one can see from Figs. (Ha)-(c), axions from sax- 
ion decay can contribute to the amount of extra radi- 
ation. However, for the considered x = 1 case, val- 
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FIG. 2. (a)-(c) Contours of AN e g = 0.26 (dash-dotted), 0.52 (dotted), and 0.78 (dashed) provided by axions from decays of 
thermal saxions and of fig/i 2 = 0.129 (solid) in the mg-TR parameter plane for gravitino LSP scenarios with m CT = mg and 
x = 1. Black (gray) curves are obtained with m\i% = 400 (500) GeV and ma = 1 (1.25) TeV. The PQ scale and the axino 
mass are set to (a) / PQ = 10 10 GeV and m- a = 2 TeV, (b) /pq = 5 x 10 10 GeV and m a = 3 TeV, and (c) / PQ = 10 11 GeV 
and ma = 6 TeV, respectively. Regions above the solid lines are disfavored by fig exceeding ficDM at the 3<r level, (d) The 
dilution factor A as a function of the reheating temperature Tr for x = 1. The black (gray) solid, dashed, and dotted lines are 
obtained with nij = 1 (1.25) TeV for the f a and ma combinations considered in panels (a), (b), and (c), respectively. 



ues of AN c g > 0.78 are almost completely disfavored 
by the f^~ P < Ocdm constraint if rrig — 1 TeV and 
m 1/2 = 400 GeV. In fact, if SUSY searches at the LHC 
point to minimum rrig and values of respectively 

1.25 TeV and 500 GeV, the rff p < ft CD M constraint 
will clearly disfavor AA^g > 0.78 and the BBN-infcrred 



posterior maxima AA^ c ff = 0.76 and 0.77 given in Ta- 
ble HI Still axions from decays of thermal saxions can 
then provide a viable explanation of AN c g < 0.52. This 
matches the 2a sensitivity of the Planck satellite mission, 
which implies that a statistically significant 3cr detection 
of AN c ff cannot be expected. 
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To explore the simultaneous viability of successful lep- 
togenesis and an explanation of, e.g., AiV of j ~ 0.52 by 
axions from decays of thermal saxions, one has to con- 
sider the minimum Tr value together with the dilution 
factors shown in Fig. [2jd) as described in (fl"2j). Indeed, 
if the minimum Tr is 10 9 GeV without the entropy pro- 
ducing axino decays, it will become almost twice as large 
in the scenarios with /pq = 10 11 GeV. Accordingly, as 
can be seen in Fig. [2jc), experimental insights on m g 
and TOx/2 will decide on such a simultaneous viability for 
1 = 1. For the lower / P q values considered in Figs.[5Ja) 
and (b), that simultaneous viability is excluded already 
with nig ~ 1 TcV and m 1/2 ~ 400 GeV. 

The described pictures changes considerably if x <C 1. 
This is shown for x = 0.2 (black) and 0.1 (gray) in 
Figs. [3[a) and (b). Here the dashed, dotted, and dash- 
dotted lines indicate AiV e ff = 0.78, 1.73, and 3.59, re- 
spectively, where the latter appears in panel (b) only. 
The solid lines show fig/i 2 = 0.129 contours. In both 
panels, m§ = 1 TeV and my 2 = 400 GeV. For both 
(a) /pq = 10 10 GeV and (b) 10 11 GeV, one finds that 
the amount of extra radiation AiV e ff that is compatible 
with the f2~ p < ^cdm constraint is now significantly 
larger than in the corresponding x = 1 cases. Not only 
the posterior maxima inferred from BBN analyses and 
AiVeff = 0.78 but also the mean inferred from precision 
cosmology, AiVeff = 1-73, are easily explained by axions 
from thermal saxions in the part of the mg-TR parameter 
plane in which < ficDM- To challenge this expla- 
nation of AiVeff = 1-73 for / PQ = 10 10 GeV and x = 0.2 
(0.1), it will require a sizable increase of the gaugino mass 
limits to m g > 1.5 TeV (1.7 TeV) and m 1/2 > 600 GeV 
(700 GeV). Without more restrictive gaugino mass lim- 
its, one sees in panel (b) that the 2a upper limit from 
precision cosmology, AiV c ff = 3.59, translates into an up- 
per limit on Tr that can be more restrictive than the one 
from f2~ p < SIcdm- 

Depending on the findings of Planck for AiV c ff, Tr lim- 
its can become considerably more severe. For example, if 
Planck does point to a favored AN c g value close to zero 
and thus to a 3a upper limit of AiV c ff < 0.78 then Tr will 
be constrained from above as indicated by the dashed line 
in both plots. In the other case in which Planck points 
to a significant amount of extra radiation, e.g., with a 
favored valued of AN c g — 1.73, the dashed lines are ex- 
pected to be close to a 3a lower limit of AiV e g > 0.78. 
These lilies will then indicate a lower limit on the re- 
heating temperature Tr . This will then point to favored 
regions in the (?ng , Tr, /pq, x) parameter space that 
will be tested further by ongoing SUSY searches at the 
LHC and resulting insights on rrig and TOi/2. 

Towards smaller x values in the range 0.1 % x < 1, 
AiV e ff increases considerably because of the later decay 
of the saxion; cf. (|20|) . At the same time, there is also 
a growing branching ratio (|22j) of the entropy producing 
saxion decays into gluons. By comparing the figii 2 = 
0.129 contours shown for mi/a = 400 GeV and /pq = 



10 10 GeV in Fig. [^a) with the ones in Fig. G^a), one 
can however see that there is only a minor additional 
dilution for x = 0.1 and 0.2 due to a — > gg decays. Also 
for /pq = 10 11 GeV, the additional dilution from saxion 
decays stays rather modest at those x values. This can 
be seen explicitly in Fig. [ljb). The additional kink on 
the solid line in Fig. ^[b) that appears for x = 0.1 at 
Tr = Tp can still be understood as a manifestation of 
this. For even smaller x values below 0.1, the dilution 
from saxion decays can become substantial, as shown in 
Fig.[TJ Together with the decreasing branching ratio (|2ip . 
this then leads to a reduction of AN c g, which can also 
be seen in (|B10I) of Appendix [Bj In fact, we find the 
maximum viable AiV e ff values for x ~ 0.1. 

For 0.1 < x <C 1, we find a much larger parameter 
space that allows for the simultaneous viability of suc- 
cessful leptogenesis working a minimum Tr ~ 10 9 GeV 
and of a sizable AiV c ff provided by axions from decays of 
thermal saxions. As can be seen in Fig. 02[a), where A 
is close to 1, T R - 10 9 GeV together with AiV c ff ~ 0.8 is 
in the allowed region when x ~ 0.1, mg ~ 30 GeV and 
/pq = 10 10 GeV. Towards larger / PQ ~'l0 n GeV, where 
significantly larger values of ma can be required, larger 
Tr and larger AiV c ff values become viable already for 
x = 0.2 in the region with mg > 30 GeV; cf. Fig. [3Jb). 

At this point, it should be stressed that contributions 
to the saxion energy density can reside in coherent os- 
cillations of the saxion field. This can give additional 
and even dominating contributions to p a and thereby to 
AiV e ff 0, GJ, H3, H HI HI]. However, these contri- 
butions depend on the initial misalignment of the saxion 
field Oj. In fact, for the considered values of m a and /pq, 
the contribution of this non-thermal source is negligible 
if c?i ~ /pQi as often assumed in the literature. 

As mentioned at the beginning of this section, we focus 
here on scenarios in which f2 ~ p provides the dominant 
part of OcDM- In Figs. [5] and [Si this holds in the region 
close to the solid lines for the respective gaugino masses. 
Moving away from those solid lines towards smaller Tr, 
there is room for additional contributions to ^cdm when 
assuming the considered gaugino masses. 

There can be a contribution to Acdm from coherent 
oscillations of the axion field after it acquires a mass due 
to instanton effects at T < 1 GeV. The resulting ax- 
ion relic density from this misalignment mechanism de- 
pends on the initial misalignment angle —n<0i<ir and 

/p Q BIMlei 

n^tf- 0.15 (^^y 76 , (34) 

where £ = 0(1) parametrizes theoretical uncertainties 
related, e.g., to details of the quark-hadron transition 
and of the T dependence of m a . Moreover, /(0f ) is the 
anharmonicity factor which satisfies f{9f) — > 1 for small 
9? — > and becomes sizable towards large 0f — > n 2 
l69j | . Here the possibility of late time entropy production 
is not included that can lead to a dilution of f^ IS /i 2 . 
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FIG. 3. Contours of AiV e g = 0.78 (dashed), 1.73 (dotted), and 3.59 (dash-dotted) provided by axions from decays of thermal 
saxions and of QqIi 2 — 0.129 (solid) in the mg-TR parameter plane for gravitino LSP scenarios with m a = mg, mi« = 400 GeV, 
and ma = 1 TeV. Black (gray) curves are obtained with x = 0.2 (0.1). The considered values of the PQ scale and the axino 
mass are respectively (a) /pq = 10 10 GeV and ma = 2 TeV and (b) /pq = 10 11 GeV and ma = 6 TeV. Regions above the solid 
lines are disfavored by a gravitino density parameter that exceeds S7cdm at the 3cr level. In the regions above and to the left 
of the dash-dotted lines, which show up only in panel (b), AN e g exceeds the 2<r upper limit from precision cosmology. 



However, already for A = 1 and Q\ ~ 1, f^ IS is only a 
minor fraction of JIcdm for /pq < 10 11 GeV considered 
in this section. In fact, with the possibility of Q\ <C 1, 
57^ can be negligible and then does not tighten the Tr 
limits indicated by the solid lines in Figs. [2] and [3l 

The contribution from decays of the LOSP into the 
gravitino LSP, 

nf SP ^ dx h 2 = m s Y LO Sps(T )h 2 /p c , (35) 

depends strongly on the LOSP type, its mass and cou- 
plings, and other details of the considered point in the 
SUSY parameter space. For the case in which the light- 
est neutralino y? i s the LOSP, the yield after freezeout 
can be sizeable |70l [7lj , 

^,o S P~d-4)xl0--( I ^ F ) ) (36) 

and can thus imply Tr constraints that are significantly 
more restrictive than those shown in Figs. [2] and [3] This 
holds even with lLOSP-diluting entropy production in 
saxion decays leading to the maximum contribution of 
A a ^sg ^ 2.5 the total dilution factor A seen in Fig.QJb). 
In contrast, for a charged slepton LOSP l\ or a sneutrino 
LOSP V\ respecting the upper limits on tolosp given in 



Table imi F LO sp is relatively small (Z^|z3, 

Tootv)' < 37) 

r„ LOS P~2xl0-"( I5 =|_), (38, 

and basically negligible already without a possible di- 
lution (fTTj) . Then the Tr limits imposed by f2~ p + 
^losp^ga' f2 CDM are very similar to the ones indi- 
cated by the solid lines in Figs. [2] and [3l 

In the considered gravitino LSP scenarios with a LOSP 
being the NLSP, the LOSP has a long lifetime before de- 
caying into the gravitino. Often such decays are found 
to take place during and after BBN. For the Xi LOSP, 
decays such as Xi ~^ Gqq can then reprocess the pri- 
mordial light elements via electromagnetic and hadronic 
energy injection. Thereby, the observationally inferred 
primordial abundances of those elements translate into 
upper limits on Ylosp that depend on the lifetime of 
the LOSP tlosp- Towards small values of tlosp which 
occur towards smaller values of ttIq, the Ylosp limits be- 
come weaker and disappear. For a Xi LOSP with mass 
m x° x — to losp an d the latter given in Table IIII1 BBN 
constraints exclude mg > 1 GeV and thereby most of 
the interesting parameter regions considered above. 

For the charged slepton and sneutrino LOSP cases, 
hadronic energy injection requires 4-body decays such 
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as l\ — > Glqq or v\ —> Gvqq and is thereby less effi- 
cient. However, a long-lived charged slepton can form 
bound states with the primordial nuclei and thereby cat- 
alyze, e.g., the primordial production of lithium-6 sub- 
stantially. This catalyzed BBN (CBBN) then imposes 
the upper limit < 5 x 10 3 s. Together with an up- 
per limit on the slepton mass mj < "t-losp ^ 300 GeV 
imposed by axino cosmology, this translates into the con- 
straint mg < 4 GeV, which again disfavors most of the 
interesting parameter regions considered above. More- 
over, the current lower limit from searches for lon g-li ved 
charged sleptons at the LHC, mj > 300 GeV [73LTt3| . is 
already in conflict with most of the mpQgp values listed in 
Table Hill In fact, the production of a long-lived charged 
slepton LOSP could leave clear signatures at the LHC 
and allow for a precise measurement of its mass. For 
example, with mj ~ 400 (500) GeV, axino-imposed con- 
straints become difficult to evade and the CBBN limit on 
tj- implies < 6 (10) GeV. Such a discovery will thus 
not be compatible with the AiV G ff explanation via decays 
of thermal saxions for Tr > 10 9 GeV. It may instead 
point to smaller T R < 10 s GeV, smaller / PQ < 10 10 GeV, 
and larger ma or to the axion CDM scenarios with an 
eV-scale axino LSP and the gravitino NLSP considered 
in the next section. 

For the sneutrino LOSP case, the presented high-TR 
explanations of additional radiation via decays of ther- 
mal saxions are still viable. The BBN constraints im- 
posed by hadronic and electromagnetic energy release be- 
come relevant only for large > 500 GeV. At smaller 



/??,- 



< m 



max 
LOSP' 



the only bound on the gravitino mass 
then results from the hierarchy mg < assumed in 

this section. In comparison to the l\ LOSP, it will be 
much more challenging to identify a sneutrino v\ as the 
LOSP and to measure its mass at the LHC. Such a mea- 
surement will allow us to test the presented scenarios in 
two ways: (i) by confonting m^ with the upper limit 
to losp imposed by the axino and (ii) by exploring the 
maximum Tr values for the maximum viable mass of the 
gravitino LSP which is then = mp r 



V. AXION CDM CASE 

In this section we consider SUSY scenarios in which 
the axion with a mass of m a ~ 6 /ieV provides the CDM 
density SIcdm^ 2 via the misalignment mechanism. The 
associated relic density D,^ is h 2 that resides in coherent 
oscillations of the axion field is given by (IMl) in the ab- 
sence of late-time entropy production. With entropy pro- 
duction after the QCD phase transition, T « 1 GeV, 
the corresponding dilution factor A has to be taken into 
account that reduces the density parameter by a factor 
of 1/A. Accordingly, n™ ls h 2 = l] CDM h 2 holds with 
/ PQ = 10 12 GeV, e.g., for /(6>f)0? - 1 and A - 1 or 
equally for /(0 2 )0? - 10 and A - 10. In fact, also in 
the situations with a sizable A ~ 30 encountered below, 
the CDM density can be explained fully by the axion 



condensate provided 6f and the associated anharmonic- 
ity factor f(9f ) are sufficiently large to compensate for 
the dilution. 

In a setting with ^cdm provided by the axion conden- 
sate, the LSP is no longer required to be a CDM particle. 
In turn, the LSP can be a very light particle such as an 
axino with ma < 37 eV, which is the scenario considered 
in this section. Such a light axino can still be produced 
thermally when Tr < Xg or decouple as a thermal relic 
when Tr > Tq. The resulting population can contribute 
to hot dark matter (HDM). In fact, the upper limit of 
ma < 37 eV is inferred from LSS constraints on HDM 
contributions in mixed models with CDM 75j. When 
relativistic, the axino population from thermal processes 

contributes a small amount of (AiV off )~ q/TP < 0.017 [zU 
to dark radiation, which is not included in our calcula- 
tions below. 

In our considerations the gravitino is the NLSP that 
is lighter than the LOSP, i.e., than the lightest sparti- 
clc in the MSSM. Thereby, the explored scenarios are 
not subject to the restrictive upper limits on Tr imposed 
by BBN constraints on hadronic/electromagnetic energy 
injection in late decays of gravitinos into MSSM parti- 
cles [7(1 [77j ■ In the R-parity converving settings consid- 
ered in this section, gravitinos can decay into axions and 
axinos only. The gravitino lifetime rg is then governed 
by the associated decay rate [zl, Wf\ 



G— >aa 



1 



(39) 



where M P = mp/y/Sn = 2.44 x 10 18 GeV is the re- 
duced Planck scale. Accordingly, gravitinos can be very 
long-lived. For example, rg ~ 10 10 s and 10 5 s for 
mg ~ 60 GeV and 6 TeV, respectively. The axions and 
axinos emitted in decays of a thermally produced grav- 
itino propulation can thereby contribute substantially to 
AAcff at late times well after BBN [1 [ll [20| . In fact, the 
time at which the smallest observable modes of the CMB 
reenter the horizon, t — 5.2 x 10 10 s, imposes an upper 
limit on Tg because of the non-observation of a significant 
release of extra radiation thereafter [l9[ . The correspond- 
ing mass limit is mg > 35 GeV. For m CT = mg, this limit 
implies that saxions decay before the onset of BBN even 
when small x values of are considered. With (f2T)]) being 
valid to a very good approximation in this section also, 
one finds r CT < 0.5 s for / PQ = 10 12 GeV and x = 0.02. 

While the mass hierarchy m a 3> ma and the La- 
grangian (|17[) now allow for the additional a — > da decay 
channel, the corresponding decay width 



PQ 



2m a 



(40) 



is suppressed by a factor of at least 32m|/m 2 r with re- 
spect to T a ^ aa given in (fT8)l and thereby negligible for 
the considered mass splittings. The saxion lifetime and 
the relevant branching ratios are thus again described 
by (HI, (HQ), and J22]), respectively. 
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As in the previous section, we encounter the two types 
of decays of non-relativistic particles. However, in the 
scenarios in the previous section, A.N e g originates from 
saxion decays only and entropy production at two very 
different times is possible. Now there are two sources 
of extra radiation, saxion decays and gravitino decays, 
which proceed at very different times, whereas entropy 
can be produced in saxion decays only. In the following 
we thus calculate 



AA cff (T) 



120 

77T 2 T,f 



Pdr(T), 



(41) 



where the energy density of dark radiation pdr includes 
contributions of axions from thermal processes in the 
early universe, of axions from decays of thermal saxions, 
and of axions and axinos from decays of thermally pro- 
duced gravitinos. As in the previous section, the possibil- 
ity of entropy production in saxion decays is taken into 
account and decays are treated beyond the sudden decay 
approximation. Thereby, we update and generalize ex- 
isting results presented in Refs. [4. [20|. For a qualitative 
understanding of our numerical results, we again refer to 
the expressions obtained in Appendix [B] 

The following Boltzmann equations describe the time 
evolution of the energy densities of gravitinos, saxions, 
and dark radiation 

pg + 3Hp e = -T dPd , (42) 
per + 3Hp a = -T a p ai (43) 
Pdr + 4Hp dT = BR(cr aa)T <J p <J + T d p d , (44) 

in the epoch well after the one in which thermal processes 
involving EWIP were efficient and when gravitinos and 
saxions from such processes are non-relativistic. Here the 
time evolution of the entropy S and the scale factor R 
are described respectively by 



1/3 



[1 — BR(cr — > aa)] T a p a 



and the Friedmann equation 
8tt 



H 1 



3ml 



(Pdr + Pa + Pq + Prad) 



(45) 
(46) 



with p ra d as given in (|29|) . 

We solve the closed set of differential equations (|42[) 
(|46p numerically. As in the previous section, we start 
at ti — 1.6 x 10 _13 s corresponding to Tj = 1 TeV with 
R[ti) = 1 GeV" 1 . However, the considered end of the 
evolution is now set to a much later time of tj = 10 12 s 
corresponding to 7/ ~ 1 eV. The initial values of the 
energy densities are given by 



p d (U) = m d Y? p s(Ti), 
p u (t i ) = m a Y^ F s(T i ), 



(47) 
(48) 
(49) 



and of the entropy by S(ti) = s(Ti)R(ti) 3 . Entropy pro- 
duction is quantified by the dilution factor A given as 
in (|33|) . Note that the contribution of the energy den- 
sity of axinos from thermal processes in the early uni- 
verse can be neglected in (|46|) at the considered times. 
Also in (|49|) . this population is neglected, which con- 
tributes at most (AA r e ff)~ q ^ TP ~ 0.017 [z^|, as mentioned 
above. The Boltzmann equation for cold dark matter 
axions and the associated contribution in (14T>|) are not 
mentioned above. In fact, including this population ex- 
plicitly leads to at most a 1-2% effect in AN c g and only 
in settings with > 10 10 s. As in the saxion treat- 
ment in the previous section, saxions and gravitinos are 
described as non-relativistic species throughout the time 
interval [ti, tf] although saxions and gravitinos, e.g., with 
m a g = 100 GeV will be relativistic at an initial tempera- 
ture of Ti = 1 TeV. This simplified treatment is justified 
since the contributions of saxions and gravitinos to the 
right-hand side of the Friedmann equation (|46[) become 
relevant only when they are non-relativistic. 

Figure BJa) presents the results of the numerical inte- 
gration for m a = m d = 100 GeV, T R = 1.3 x 10 10 GeV, 
/pq = 10 12 GeV, and universal gaugino masses at the 
GUT scale of m 1 / 2 — 400 GeV, which is compatible 
with mg = 1 TeV at collider energies. The time evo- 
lution of R?p is shown for saxions (dashed), gravitinos 
(dash-dotted), dark radiation (dotted), and other radi- 
ation (solid), where black and gray lines refer to x = 1 
and 0.02, respectively. The evolution of entropy S is not 
shown. For x = 1, it is simply a horizontal line and 
A = 1. In the case with x = 0.02, it shows an increase 
by a factor of A ~ 3 when the saxion decay occurs. The 
latter dilution factor can be inferred also from the differ- 
ence of the two solid curves. Here one can see that the 
energy density of the universe can be dominated by non- 
relativistic saxions just before/during their decay, which 
indicates an early intermediate matter-dominated epoch. 
In fact, such an epoch can be even more pronounced to- 
wards larger Tr and/or smaller m a and thereby lead to 
significantly larger A values, as illustrated in Fig. Sfb). 

In Fig. HJb) the Tr dependence of the dilution factor 
A is shown for x = 1, 0.2, 0.1, and 0.02 by the solid, 
dashed, dotted, and dash-dotted curves, respectively. 
Black (gray) lines refer to m a = ra a = 50 (100) GeV, 
whereas all other parameters are as in panel (a). The 
Tr dependence results from the one of v ff cq/TP , which 
explains the kinks at = encountered already in 
the previous section. Again there is an increase of A to- 
wards small x due to larger values of t g and BR(cr — > gg) . 
For x < 0.1 and towards large Tr, > Tq, A can now be 
much larger than in the previous section because here 
/pq = 10 12 GeV. The latter implies a larger saxion life- 
time ([20]) . whereas Y£ q is independent of /pq; see also 
approximations (|B6|) and (|B7[) in Appendix |B| 

Let us now turn to AN c g. The dotted line in Fig. 0Ja) 
illustrates that there are two substantial contributions 
at very different times, as advertised above: AJV"^"" 
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FIG. 4. (a) Time evolution of the energy per comoving volume, R 3 p, of saxions (dashed), gravitinos (dash-dotted), dark 
radiation in the form of axions and axinos (dashed) and other radiation (solid). Here m CT = mg = 110 GeV, ma < 37 eV, 
mi/2 = 400 GeV, T R = 1.7 x 10 10 GeV, and / PQ = 10 12 GeV. The initial value of the scale factor is set to Ri = 1 GeV" 1 
at the temperature Ti = 1 TeV corresponding to a time of ti — 1.6 x 10~ 13 s. Black (gray) lines refer to the case with x — 1 
(0.02). (b) The dilution factor A as a function of the reheating temperature Tr for x = 1, 0.2, 0.1, and 0.02 shown by the 
solid, dashed, dotted, and dash-dotted lines, respectively, in axino LSP scenarios with the gravitino NLSP. Black (gray) lines 
are obtained with m CT = 50 (100) GeV, whereas all other parameter are as in panel (a). 



residing in axions from decays of thermal saxions and 
AN^f^ aa residing in axions and axions from decays of 
thermally produced gravitinos. For x = 1, the param- 
eters are chosen such that the early contribution agrees 
basically with the posterior maxima of the BBN analyses 
quoted in Sect. HI i.e., AN°£> aa ~ 0.75. Moreover, this 
parameter choice comes with an additional late contribu- 
tion of AN^f*™ = 1.01 so that the sum of both agrees 
basically with the mean from the precision cosmology 
study quoted in TableU i.e., AN^ aa +AN^ a " = 1.76. 
Indeed, although not statistically significant, we find it 
tempting to suggest the difference between that mean 
and the mentioned posterior maxima as a first indica- 
tion towards the realization of such a setting in na- 
ture [4j. Further supporting features are the axion con- 
densate and thermal leptogenesis providing natural ex- 
planations of CDM for /pq = 10 12 GeV and of the baryon 
asymmetry for Tr > 10 9 GeV, as already emphasized in 
Refs. USEi. 

In the following we systematically explore AN e s con- 
tributions in settings with /pq = 10 12 GeV and large Tr. 
In addition to the latter two features mentioned above, 
the saxion energy density residing in coherent saxion os- 
cillations with <Ji ~ /pq is negligible with respect to the 
one from thermal processes in that parameter region [l8j . 
Our results are presented in Figs. [5] and [5] in the mg-Tp 
parameter plane for m a = rn^, < 37 eV, and univer- 



sal gaugino masses at the GUT scale of mx/2 = 400 GeV 
or 600 GeV. The latter is compatible with rrig = 1.5 TeV 
at collider energies. The region with rg < 5.2 x 10 10 s 
is not considered and indicated by a vertical gray dotted 
line at mg ~ 35 GeV. 

In Fig.[SJa) the solid black (gray) lines show AN° ff +aa + 
AjSfG^aa = q 7s and 3 59; as labeled, and the dashed 

lines ANg~^ ad = 3.59 for my 2 = 400 (600) GeV and 
x = 1. To allow for a comparison, the diagonal dotted 
line indicates AN^f > ' aa = 3.59 as obtained for mg = 
1 TeV with the existing result of Ref. [2(| based on the 
sudden decay approximation. The difference with respect 
to the corresponding dashed line is due to the sudden de- 
cay approximation, which overestimates AN c g by about 
13%, and the omissions of electroweak and spin-3/2 con- 
tributions in the gravitino yield Y~ p used in Ref. [2fl . 
Including the electroweak contributions increases Y~ p 
by about 20% at mg ~ 35 GeV, while the importance 
of the spin-3/2 components becomes much more pro- 
nounced towards larger mg. Comparing the respective 
dashed and solid lines, we find that AN^ aa contribu- 
tions lead to an additional sizable AN c g increase. In fact, 
for m a > 100 GeV, they tighten the upper limit on 
imposed by the 2a upper limit from precision cosmology, 
AN e s < 3.59, by up to almost one order of magnitude. 
The kinks on the solid contours at Tr = Tg reflect again 
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FIG. 5. Contours of A7V c ff provided by axions from decays of thermal saxions, AN^ aa , and by axions and axinos from decays of 
thermally produced gravitinos, ANtff^ aa , in the wi^-Tr parameter plane in axino LSP scenarios with the gravitino NLSP, where 
m CT = mg, ma < 37 eV, and /pq = 10 12 GeV. In panel (a), x = 1 and black (gray) contours refer to mi/2 = 400 (600) GeV. 
Here we show solid contours of AN^ aa + AN^ a& = 0.78 and 3.59 and dashed contours of AN^ ali = 3.59. The diagonal 
dotted line indicates the latter as well but as obtained with the AN§^ aa estimate from Ref. [2Q|. In panel (b) we show 
solid, dashed, dotted, and dash-dotted contours of AN^ aa + AN§^ ai = 3.59 for x = I, 0.2, 0.1, and 0.02, respectively, and 
m i/2 = 400 GeV. The regions above these contours are disfavored by precision cosmology at the 2a level. The vertical dotted 
line indicates the lower limit on mg from the requirement rg < 5.2 x 10 10 s in both panels. 



the Tr dependences of the saxion yield. To see this more 
explicitly and for further comparison, we refer to Fig. 6 in 
Ref. [H|, where AA^^ ao only is presented as obtained in 
the sudden decay approximation. Also (|B10|) and (|B13p 
in Appendix [B] of this work are approximate analytical 
expressions respectively for AN^* aa and AN^~ > ' aa that 
are based on the sudden decay approximation. 

Figure E]Ja) demonstrates how the AN e ff contours will 
move if LHC experiments point to rrig > 1.5 TeV and 
thereby to TO1/2 > 600 GeV. These changes are gov- 
erned fully by AN^ ad whereas AN^ aa j s not affected. 
At this point, we should stress that a collider measure- 
ment of the LOSP mass tolosp will limit from above. 
While the chosen mi/ 2 values can imply an mLOSP value 
that is well below 1 TeV, we refrain from presenting such 
an upper limit for mg since it will depend strongly on 
other details of an assumed SUSY model as well. 

The difference between the solid AA c ff = 3.59 and 
AA e ff = 0.78 contours illustrates the possible impact 
of results from the Planck satellite mission. If Planck 
should not find any statistically significant hints for ex- 
tra radiation, then the contour AA ff = 0.78 can provide 
the new upper limit on Tr at the 3<7 level. For x = 1, the 
viability of Tr > 10 9 GeV will then depend on tolosp 
and on other LOSP-related cosmological constraints dis- 



cussed below. 

Let us now turn to the case of x < 1. Figure[5fb) shows 

AJV^"" + AN^ ah = 3.59 contours for x = 1 (solid), 
0.2 (dashed), 0.1 (dotted), and 0.02 (dash-dotted), where 
m i/2 — 400 GeV. Corresponding dilution factors A have 
already been shown in Fig.2£b) and discussed thereafter. 
The x dependence of AN^f^ aa results fully from the one 
of A so that this contribution decreases towards x —> 0. 
In contrast, for the same reasons as in the previous sec- 
tion, AJV"^"" increases towards smaller x in the interval 
0.1 < x < 1, reaches its maximum at x ~ 0.1, and de- 
creases thereafter, i.e., towards smaller x < 0.1. The 
latter behavior transfers to AN^ aa + AJV^"", as can 
be seen in Fig. EJb). Here the most restrictive upper Tr 
limit is found for x = 0.1 and the most relaxed one for 
x = 0.02, where AN^ aa + AN^ ad ~ AN^ a ". 

In Fig. E] we explore how AN^ aa + AAj^ aa (solid) 

can emerge as a composition of a late AN^f^ aa (dashed) 
and an early AN^ aa (dotted) for (a) x = 1, (b) 0.2, 
(c) 0.1 and (d) 0.02. In each panel, we consider m-i/2 = 
400 GeV and show gray and black contours of AN e g = 
0.78 and 1.73, respectively. For the reasons discussed 
above, i.e., taken the posterior maxima and the mean 
listed in Table U at face value, we are particularly in- 
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FIG. 6. Contours of AN a g provided by axions from decays of thermal saxions, AN^ aa (dotted), by axions and axinos from 
decays of thermally produced gravitinos, AN^f^ aa (dashed), and by the sum of both, AN^ aa + ANg^ aa (solid), in the 
tUq-Tb, parameter plane in axino LSP scenarios with the gravitino NLSP, m CT = rreg, ma < 37 eV, /pq = 10 12 GeV, and 
mi/ 2 = 400 GeV. Black (gray) lines indicate AN c f{ — 1.73 (0.78). In each of the four panels, a different x value is considered: 
(a) x = 1, (b) 0.2, (c) 0.1, and (d) 0.02. The vertical dotted line is as in Fig. [5] 



trigued by an early AN c g ~ 0.78 from saxion decays and 
a total late AN e g ~ 1.73 after gravitino decays. The pa- 
rameter points that allow for such a composition are the 
ones at which the gray dotted and the black solid lines 
intersect. Accordingly, this composition is possible in all 
four panels, i.e., for x = 1, 0.2, 0.1, and 0.02. Note that 
there is no AN^* aa = 1.73 contour in panel (d). In fact, 
the largest possible value of that contribution in the pa- 



rameter region shown in that panel is AN^ aa ~ 0.93 
which occurs for m a ~ 400 GeV and T R > Tg. 

In light of the BBN uncertainties with respect to an 
early AN c g, it should be emphasized that compositions 
different from the one considered in the preceding para- 
graph are possible as well, as discussed at the end of 
Sect. (Til It is insightful to consider the contours in Fig. [S] 
also in light of the different possibilities of the outcome 
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of the Planck satellite mission outlined in Sect. [TTJ For 
example, if the new mean of a late AiV e g- turns out to 
be close to 0.78, the AN c g = 1.73 contour would then be 
close to a 3cr upper limit on Tr, whereas combinations 
near the gray solid lines would be the favored ones. On 
the other hand, if that mean is close to 1.73, then the 
AiVeff = 0.78 contour will be close to a 3<r lower limit on 
Tr, whereas combinations near the black solid lines will 
be the favored ones. In fact, we consider it best to revisit 
the explored settings once the relevant Planck results are 
published. 

In an assessment of the simultaneous viability of ther- 
mal leptogenesis and a certain AiV e g- composition, the 
corresponding dilution factor A has to be taken into ac- 
count in the same way as in the previous section. While 
this factor can now be much larger, the current upper Tr 
limit imposed by AiV e ff < 3.59 still allows for that simul- 
taneous viability even for x = 0.1 when > 50 GeV; 
cf. Figs. EJb) and EJb). Planck results may lead to a 
more restricitve AiV c ff limit. Nevertheless, even with 
AiVeff < 0.78, thermal leptogenesis can remain viable 
for x = 1 and also for smaller x provided can be suf- 
ficiently large. Because of the assumed hierarchy in this 
section, mg < mLOSP, a measurement of ttilosp can 
thus challenge that simultaneous viability, in particular, 
for x < 0.1 and larger mi/2- 

In the considered situation with the axino LSP and 
the gravitino NLSP, the LOSP is again a long-lived par- 
ticle. However, in contrast to the gravitino LSP set- 
ting in Sect. IIV1 it can not only decay into gravitinos, 
LOSP — > GX, but also into axinos, LOSP — > aX, where 
the relative importance is governed by mg and /pq. For 
/pq = 10 12 GeV and mg > 35 GeV, the decay into the 
axino is the dominating one, i.e., T(LOSP — > aX) 
r(LOSP -> GX). Thereby, the LOSP lifetime can be 
significantly shorter than in the previous section so that 
the (C)BBN constraints related to a late decaying LOSP 
described at the end of that section can be evaded. In 
fact, the charged slepton LOSP is now a viable possibil- 
ity, which is particularly attractive since it could appear 
as a quasistable charged massive particle in collider ex- 
periments. For example, if the LOSP is the lightest stau 
with mf x > 300 GeV, there is indeed no limit on the grav- 
itino mass other than mg < for /pq <5x 10 12 GeV 
and already with A = 1 [75[. Late time entropy pro- 
duction in saxion decays with x -C 1 can dilute Ylosp 
as described by (jlip with a sizable A > 1 and thereby 
imply even more relaxed constraints. The bino-like neu- 
tralino LOSP situation was considered in Ref. [78[ and 
found to be viable for /pq ~ 10 12 GeV as well. This work 
accounted for entropy production in saxion decays also 
but did not address the production of extra radiation. 
Similarly, the sneutrino LOSP situation is expected to 
be viable in the considered settings. 

Very light axinos with m a < 37 eV emitted in LOSP 
decays contribute only negligible amounts to the density 
parameter and to extra radiation AiV e ff. This can be seen 



by evaluating 

n L a OSP ^~ aX h 2 ~ y/pl i0 +mlY L osps(T Q )h 2 /p c , (50) 

with the present momentum of these axinos Pa,o as ob- 
tained in the sudden decay approximation. For the 
AN^ SP ^ aX contribution, this is shown explicitly for 
the stau LOSP case in Sect. 4.3 of Ref. [7j|, which can 
easily be generalized to other LOSP candidates. Axions 
and axinos from gravitino decays are still relativistic to- 
day for the considered values of m a and mg. Accord- 
ingly, their contribution to the density parameter can be 
expressed in terms of AJV^"": 

n G^ad h 2 + QG^aa h 2 = 57 x 1Q -6 & N G^aa _ ^ 

This clarifies that the shown constraints will neither be 
tightened by AiV L ff OSP ^ 5X nor by contibutions to the 
density parameter described by (|5T)|) and (fSTj) . 

The density parameter of axinos from thermal pro- 
cesses in the early universe is given by 

^T /TP h 2 - V /(^ h ) 2 +- 2 1T /TP s(T )h 2 /p c , (52) 

where the average momentum of thermal axinos today is 
given 

by (Pto) = 3.151 Ta,o with the present axion tem- 
perature of T aj0 = [3* s (7o)/228.75] 1 / 3 T ~ 0.06 meV. 
Expression ((52)) relies on the fact that not only thermal 
relic but also thermally produced axinos show basically 
a thermal spectrum. When comparing (p~ h ) with the 
axino mass, one finds that this axino population can still 
be relativistic today but only when m a < 0.2 meV. A 
similar comparison shows that such a thermal axino pop- 
ulation with, e.g., m a < 4 eV was still relativistic at t = 
5.2 x 10 10 s. As mentioned at the beginning of this sec- 
tion, this axino popoulation is dark radiation when rela- 
tivistic. When non-relativistic, the energy density of this 
axino population and its contribution to the density pa- 
rameter are governed by m a . This allows for quantifying 
the HDM constraint as m a < 37 eV [75j. Interestingly, 
the considered axion CDM hypothesis will continue to 
be probed by the direct axion search experiment ADMX 
exactly in the region around /pq ~ 10 12 GeV [79j]. A 
discovery of axions in this search together with a statis- 
tically significant evidence for extra radiation from the 
Planck satellite could therefore point towards the real- 
ization of one of the settings considered in this section. 
Further support in favor of those settings (and against 
the ones considered in Sect. lIV[) would be the discovery of 
a long-lived charged slepton LOSP at the LHC. It could 
even be that future cosmological analyses find hints on 
the time of the release of extra radiation. Such a release 
may manifest itself in the perturbation spectrum so that 
precision cosmology might help to assess the lifetime of 
the gravitino whose late decays produce dark radiation 
at times before 5.2 x 10 10 s. Another strong hint for the 
scenarios considered here would be the confirmation of 
extra radiation prior to BBN and a significant difference 
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between that amount with respect to the one at much 
later times. This would require advances in BBN-related 
studies. In particular, this calls for new high quality spec- 
tra from extragalactic HII regions that should allow for 
a significantly more precise determination of AN c g prior 
to BBN 



VI. CONCLUSION 

We have explored two scenarios of hadronic axion mod- 
els [H| in R-parity conserving SUSY settings: (i) a 
gravitino LSP scenario with a heavy axino at the TeV 
scale, and (ii) a scenario with a light axino LSP at the eV 
scale and the gravitino NLSP. Both scenarios are found to 
allow for consistent explanations of extra radiation and 
CDM and for a high reheating temperature Tr of up to 
about 10 9 GeV or 10 11 GeV, respectively. Testable cases 
have been outlined that may still allow for the high Tr 
values required by successful thermal leptogenesis with 
hierachical heavy Majorana neutrinos |29j . 

In the gravitino LSP scenario, CDM resides in the 
gravitino and AN c g is explained by thermal saxions 
which decay into axion pairs prior to BBN. We have 
shown that up to AN e s — 0.8 can arise naturally for 
10 10 GeV < f PQ < 10" GeV and 10 7 GeV < f R < 
10 9 GeV. This finding requires that the gluino mass m§ 
is close to the current experimental limit of about 1 TeV. 
For a larger rrig = 1.25 TeV, we have demonstrated that 
smaller values of AN e g ~ 0.5 remain viable. Viability of 
larger AN c g (i.e., above 0.8 or 0.5) is found to require 
a more suppressed saxion- axion coupling, i« 1, with a 
maximum AN c g occuring for x ~ 0.1. There a statis- 
tically significant AN e g > finding at Planck can even 
lead to lower limits on Tr and thereby point to favored 
regions in the parameter space. 

For compatibility of the presented gravitino LSP case 
with cosmological constraints, the axino must be heavy, 
rria > 2 TeV, so that it decays prior to the decoupling 
of the LOSP from the thermal bath. For the high Tr 
values considered, such a heavy axino can still be pro- 
duced very efficiently in thermal processes in the early 
Universe. Primordial axinos can thereby contribute sig- 
nificantly to the total energy density just before decay- 
ing into gluinos and gluons. Calculating the associated 
entropy production, we obtain dilution factors of up to 
fca-*99 r^, 2 that affect the abundances of gravitinos, sax- 
ions, and axions produced in thermal processes well be- 
fore axinos dominate the energy density. Since also a 
baryon asymmetry generated prior to that epoch is di- 
luted by the same factor, about twice of the observed 
value is needed prior to that dilution. Within the frame- 
work of thermal leptogenesis, this implies that the usually 
required Tr ~ 10 9 GeV j2!j now has to be basically twice 
as large [33, H3- For x = 1, we find this to be viable for 
/ PQ = 10 11 GeV, m~ g ~ 1 TeV, and m~ a > 6 TeV when 
AN e g < 0.5. Towards small x <C 1, the decay of thermal 
saxions into gluons can lead to an additional sizable di- 



lution factor of up to A a ^* 99 ~ 3. This can dilute even 
the yield of the LOSP after decoupling from the ther- 
mal plasma and prior to decay and thereby weaken BBN 
contraints related to late decaying LOSP. Moreover, for 
0.1 < x -C 1, we have found that a significant part of the 
parameter space will allow for the simulaneous viability 
of thermal leptogenesis, a sizable AjV e ff provided by ax- 
ions from decays of thermal saxions, and SIcdm residing 
thermally produced gravitinos. 

In the scenario with the light axino LSP and the 
gravitino NLSP, CDM resides in axions from the mis- 
alignment mechanism, which provides naturally ri^ IS — 
Qcdm for fpQ ~ 10 12 GeV. Remarkably, the ongoing di- 
rect axion CDM search by ADMX [79| is sensitive in ex- 
actly that /pq range and may find signals supporting this 
CDM explanation in the near future. We have demon- 
strated that there are now two sources for a possibly sub- 
stantial AN c g that work at very different times: thermal 
saxions that decay into axion pairs prior to BBN and 
thermally produced gravitinos that decay into axions and 
axinos well after BBN and before 5.2 x 10 10 s. Accord- 
ingly, we find different possibilities to explain possible fu- 
ture AiVeff results of the Planck satellite mission within 
this scenario. If Planck finds a statistically significant ex- 
cess of dark radiation, e.g., of AN e s ~ 1.7±0.8 at the 3a 
level, a natural explanation will be the composition with 
an early AN^ aa ~ 0.8 residing in axions from saxion 

decays and an additional late AN^f^ aa ~ 0.9 residing in 
axions and axinos from gravitino decays. In fact, this hy- 
pothesis is supported by posterior maxima found in BBN 
likelihood analyses [3j, |4[ based on the recent Y p studies 
of Refs. [HQ. However, without more precise BBN lim- 
its for AiVoff, which may indeed be difficult to obtain in 
light of the systematic uncertainties Q , there remains a 
significant uncertainty with respect to the amount of an 
early AN e ff. Accordingly, a possible Planck finding of a 
late AN g ~ 1.7 can result equally well either dominantly 
from late G — > aa decays for <C 100 GeV or domi- 
nantly from a — y aa prior to BBN for m^, ^> 100 GeV. 

Our refinements with respect to Refs. 3,H3] have been 
found to have the following effects. By treating decays 
beyond the sudden-decay approximation, the resulting 
AiV e )f values decrease by about 10%. Moreover, with 
the gravitino yield that accounts for the gravitino-spin- 
3/2 components and for electroweak processes, previously 
neglected contributions to AN^[^ aa are included which 
become sizable for > 100 GeV. Together with the 
contributions from a — )■ aa decays, this results in signif- 
icantly larger AN e R values in that region. In turn, our 
new upper bounds on Tr for = O(100 GeV) are sub- 
stantially more restrictive than previously expected for 
x = 1 in the axino LSP case with the gravitino NLSP. 
Towards a more suppressed saxion-axion coupling with 
x ~ 0.1, we find even larger AN c g values that further 
tighten the Tr limits significantly. Even smaller values 
of x <C 0.1 have been found to come with significant di- 
lution factors of up to A a ^* g9 ~ 30. Those can reduce 
AiVeff and in turn relax the upper bounds on Tr but 
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have to be included in an assessment of the viability of 
thermal leptogenesis. 

We have discussed ways in which the different expla- 
nations of a potentially sizable AN e g will be narrowed by 
ongoing SUSY searches at the LHC. Particularly impor- 
tant will be new limits on nig or measurements thereof. 
The gluino mass governs the thermally produced grav- 
itino yield and thereby limits that relate to this quan- 
tity. Upper limits on Tr thus become more restrictive 
for larger nig. In fact, for rrig > 1.1 TeV, we find that 
a AN e ff > 0.8 explanation by axions from thermal sax- 
ions becomes incompatible with the then too restrictive 
Qq < ficDM constraint in the gravitino LSP case. In the 
alternative axino LSP case, on the other hand, it is the 
AiV e g-imposed limit that becomes more restrictive to- 
wards large m g in the small mg region where the decay 

G — > aa contributes significantly to AiV e gf. 

Other relevant LHC findings will be a discovery of 
the lightest sparticle within the MSSM (i.e., the LOSP), 
its identification, and a measurement of its mass. In 
both of the considered cases, the LOSP mass limits mg 
from above. Moreover, the LOSP is expected to be 
long-lived so that additional restrictive cosmological con- 
straints can occur depending on the nature of the LOSP. 
For example, for a long-lived charged slepton LOSP, 
which has to be heavier than about 300 GeV W^ . 
CBBN constraints can disfavor the presented gravitino 
LSP case [8l| . Remarkably, such an LOSP is found to be 
compatible with the light axino LSP scenario with the 
gravitino NLSP The discovery of such an LOSP 

could thus become an important additional hint in fa- 
vor of the latter scenario. In the gravitino LSP scenario, 
BBN constraints associated with hadronic energy injec- 
tion disfavor the possibilities of a neutralino LOSP or 
a colored LOSP as well. Nevertheless, that scenario is 
found to be viable with a sneutrino LOSP. One will then 
face the challenge to identify a long-lived sneutrino as 
the LOSP |82h85I |. which will be a much more difficult 
task than the identification of a long-lived charged slep- 
ton LOSP. 

In summary, we find the presented scenarios appeal- 
ing from the cosmological point of view and intriguing 
with respect to their testability. In light of those fea- 
tures, it will be interesting to see ways in which model 
building can allow for the suggested mass spectra and the 
large splittings between the axino mass and the masses 
of the saxion and the gravitino. With upcoming new 
results from the direct axion dark matter search experi- 
ment AD MX, the Planck satellite mission, and the LHC, 
it will be exciting to see further hints for or against the 
viability of the considered scenarios soon. 

Appendix A: Thermal Axino Production 

Let us present some of the details of the calcula- 
tion that lead to our update of thermally produced ax- 
ino yield ([7]). To obtain a finite result in a gauge- 



Process H: qt + g a —¥ q~j + a 

m „ , % Qi . , 'h 




Process J: qi + q~j — > g a + a (Crossing of H) 

FIG. 7. The 2 — > 2 processes of axino production affected by 
including the quartic axino-squark-antisquark-gluino vertex 
described by the second term in the third line of @. Process 
H is also possible with antisquarks, replacing qij by qij. 



invarant treatment, we rely on systematic field theoret- 
ical methods such as HTL resummation [4l[ and the 
Braaten-Yuan prescription [42j exactly as applied in 
Ref. [iol ]. However, we now include the quartic axino- 
squark-antisquark-gluino interaction described by the 
second term in the third line of (J5J, which was not con- 
sidered in Ref. as pointed out in Ref. 'A'.) . 

Following [4Cj and the methods referred to therein 
closely, we split the thermal production rate into a soft 
part, which involves soft gluons with momentum trans- 
fer of order g s T, and a hard part, in which no soft 
gluon exchanges occur. The soft part is not affected 
by the additional vertex. In the hard part, this vertex 
contributes additional Feynman diagrams in the process 
Qi + 9 a —> q~j +ci and its crossing qi + <jj — > g a +a labeled 
respectively as H and J in |4(J . Figure [7] shows the com- 
pleted set of Feynman diagrams for process H and lists 
process J as its crossing. The corresponding squared ma- 
trix elements \Mi\ 2 for a single chirality and with sums 
over initial and final spins read 

|MH|2/ l2j^ = - 2 ( t + 2s + 2 T)'^ 2 ' < A1 > 

|Mj|2/ T28^- 2 ( S + 2< + 4) |T - 12 (A2) 

and replace the respective entries in Table 1 of Ref. ■ 
Other entries in that table are not affected. The Man- 
delstam variables are given by s = (P\ + P2) 2 and 
t = (Pi — -P3) 2 , where the particle four-momenta Pi refer 
to the particles in the order in which they are written 
down above and in Fig. [7] 

Grouping the processes into different classes (depend- 
ing on the number of external bosons and fermions) and 
weighting the matrix elements with their respective mul- 
tiplicities and statistical factors, we find that the sums of 
the corresponding squared matrix elements |Mbbf| 2 and 
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|M BF b| 2 given in (3.7) and (3.8) of Ref. 40] change to: 

I^bbfI 2 / 2£ 2 \ , N / a \ 

1 s + 2t+ — )(N c + n f ) + 4snf, (A3) 



9^-1) 



|M B fb| 2 



"32^/^" 



-t-2s 



2£ 
t 



{N c + nf)- ttrif, (A4) 



where the notation of the above reference is adopted with 
N c = 3 denoting the number of colors and nj = 6 the 
number of color triplet and anti-triplet chiral multiplets. 
Further simplifications lead to 

IMbbfI 2 9 i[ N *7, 1] [\M 3 \ 2 (N c + n f ) - 2\M 2 \ 2 n f ] , 



|.U„™l a = 1- c . .„ 1 [mi 2 (N c + n f ) - 2\M 2 \ 2 n f 



32^/r 2 



PQ 



where |Mi 



-t - 2s - (2s 2 /t) 1 \M 2 f 



(A5) 

•/]> 
(A6) 
t, and 



|M 3 | 2 = t 2 /s. Using this instead of Eqs. (3.12) and 
(3.13) of Ref. [40], our resut for the hard part of the 
thermal production rate shows a prefactor of — 2n/ in- 
stead of — 3n//2 in the fourth line of (E.l) in [|o| but 
otherwise agrees with that equation. 

By adding the soft and hard parts of the thermal pro- 
duction rate and by integrating the resulting total ther- 
mal production rate over the energy of the produced ax- 
ino, we arrive at the collision term 



W h (T) = 



(N 2 - 1) 3C(3) ffs 6 T 6 



ft 



PQ 



In 



4096tt 7 



(A7) 



1.647 T z 



(N c + n f ) + 0.5781 n, 



with m 2 = g 2 T 2 (N c + n/)/6 denoting the squared SUSY 
thermal gluon mass. The collision term enters the Boltz- 
mann equation, ha+3Hria — Wa, that describes the time 
evolution of the axino number density. Integrating this 
equation as described in [4(|, we get for the thermally 
produced axino yield 



Wg(T R ) 

s(T R )H(T R ) 



(A8) 



and thereby expression ([7]) given in Sect. Hill In summary, 
the constant in the logarithm in the expression for Y^ p 
changes from 1.211 in (E.3) of ^ to 1.271 in Q when 
including the quartic axino-squark-antisquark-gluino ver- 
tex. Accordingly, in the i?-parity conserving axino LSP 
scenarios considered in Ref. poj . the density parameter 
of thermally produced CDM axinos changes to 



n d h 2 = 5.5 5 °ln 



1.271 

9s 



0.1 GeV 



AO 11 GeV V 

^ / PQ J ^10 4 GeV 



(A9) 



Nevertheless, the qualitative statements and plots of 
Ref. [40| are only mildly affected by this correction. 



Appendix B: Approximations for A and AiV e fr 

Here we provide expressions that describe approxi- 
mately the numerical results obtained in Sects. ITVl andlVl 
The presented expressions help to understand the qual- 
itative behavior of those results and their dependencies 
on quantities such as /pq, x, Tr, to ct , mg, mg, and rrig. 

We start with the dilution factor A based on the cor- 
responding considerations in Ref. [86[ . The equation de- 
scribing the change in entropy due to the decay of a single 
non-relativistic species ip into relativistic particles that 
rapidly thermalize reads 



- 27r 2 \V3 
S ^s = R 4 {—g*s) T^, 



(Bl) 



which is the basis for (|2"T|) and (|4"5j) in the main text. Here 
PtP and are the energy density and the total decay 
width of ip, respectively, and ip is assumed to decay fully 
into rapidly thermalizing particles. By integrating (|B1[) . 
one arrives at [8a 



If = . 

Si 



dt' 



27T 2 

45 



■g*s 



1/3 



R(t>) 



-Ttb t 



. , (B2) 

where the subscript i refers to the respective quantities at 
the initial time ti. This time U can differ from the value 
used in our numerical calculations in Sects. ITVl and IVl In 
fact, the main contribution to the integral comes from the 
time interval around = 1/r^, so that, e.g., ti — O.Olr^ 
is sufficiently early to obtain a good precision. 

To solve the integral in (|B2|) . one needs to know the 
evolution of the scale factor. This is described by the 
Friedmann equation and therefore depends on the en- 
ergy content of the Universe at the relevant times. For 
the following two limiting cases, an approximate solution 
for A, the ratio of the entropy before and after the decay, 
can be obtained analytically. 

When the energy density of the non-relativistic par- 
ticle ip prior to its decay dominates the one of the Uni- 
verse, matter dominates so that R cx t 2 / z and [8a ] 



Marge 



1.83 (g* 



\l/4. 



(r^mp)V2' 



(B3) 



where (g*s) denotes a suitably averaged value of g^s over 
the integration interval. If g*s does not change signifi- 
cantly around t ~ t^, (g*s) — <?*s( T v>) gives a reasonable 
approximation. The subscript "large" in (|B3I) is used be- 
cause of the large dilution factor, A 3> 1, encountered 
in such situations and to indicate the correspondingly 
limited applicability range of (|B3[) . 

When the energy density in radiation dominates the 
one of the Universe prior and during the epoch in which 
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tp decays, R oc t 1 ! 2 and [If 



A sman ^ 1 + 1.61 



<5*5) 1/3 



5*s(^) 1/12 (r^m P )V2 



(B4) 



Again one obtains a good approximation with (g*s) — 
9*s{jip) if <?*S is (basically) constant in the relevant in- 
terval. The subscript "small" in (|B4[) indicates that its 
applicability is limited to settings in which A is not much 
larger than one. 

Let us turn to the case of entropy release from axino 
decay considered in Sect. IIV1 For all of the parameter 
points examined in this work, p^ < p 13 ,d and thus the 
corresponding dilution factor is A a ^ 99 = 0(1). Con- 
sequently, we can use (|B4j) to approximate the dilution 
factor from entropy release in axino decays. Indeed, for 
axinos from thermal processes with BR(a — > gg) ~ 1, 
our numerical results - shown e.g. by the solid line in 
Fig. [lib) - are well approximated by 



A™^l + 2.3xKr 2 



2 TeV\ 1/2 
ma J 



/pq 




-3/2 



1 - 



Y, 



iq/TP N 



io- 



10 10 GeV, 

g*s(ra) 1/3 
5 »s(0.01r a )Vi2' 

(B5) 



where ti = O.Olr^ is used as suggested above. 

Saxions can decay both into inert radiation and into 
rclativistic particles that rapidly thermalize with the re- 
spective branching ratios (f2"Tj) and (|22|) governed by x. 
For x > 0.1, A a ^ 99 = 0(1). Accordingly, after account- 
ing for BR(er — » gg), (|B4|) can be used to approximate 
the dilution factor due to entropy release in decays of 
saxions from thermal processes. Our numerical results - 
shown e.g. by the dashed and dotted lines in Fig. 0Kb) - 
are indeed well described by 



1 + 1-03x10 



_ 2 / 100 GcV \ 1/2 / / PQ 
\ m a ) ^10 10 GeV 

^° q/TP \ <Mt ct )V3 



(a 2 + 0.5x 2 7r 2 ) 3 / 2 I 10- 3 /5,s(0.01t (T ) 1 /i2' 



(B6) 



where U = 0.01r ff . For x = 0.02, / PQ = 10 12 GeV, and 
Tr > 5 x 10 10 GeV in the axion CDM scenario, A a ^ 99 > 
10 is possible and there best described by using (|B3|) . 
Setting x — in Y a and BR(cr — > gg), we then obtain 



large 



100GcVy /2 




/PQ 



10 12 GcV 



g*s{r a ) 

10.75 



1/4 



(B7) 



Fig.SJb) - in that parameter region with large A. For an 
approximate treatment of entropy production in saxion 
decays, see also Ref. |37j |. 

In settings in which two non-relativistic species decay 
at different times and thereby produce entropy at dif- 
ferent times, the total dilution factor A is given by the 
product of the individual dilution factors Aj. This oc- 
curs, e.g., for x < 1 in Sect.QVl where A = A d ^ 99 A a ^ 99 . 
There the product of (|B8[) and (|B10|) describes approxi- 
mately the dashed and dotted curves in Fig.[TJb). 

To arrive at approximate expressions for AN e g, we 
work with the sudden decay approximation as, e.g., in 
Refs. Hill. The decays that can lead to extra radiation 
are thus approximated to proceed exactly when cosmic 
time is equal to the lifetime of the decaying species. The 
contribution to AN e f{ of axions from decays of saxions 
from thermal processes is then given by 



AN^ aa 



(T) 



120 



9*s(T) 



4/3 



T 



7ir 2 T* [g*s{T a 
x BR(a -> aa)p^' TV {T a )/A 

with a temperature at the decay time t — r a of 

1/2 



(B8) 



T„ ~ 10.6 MeV 



10 10 GeV \ 
/pq / 



2a| 

7T 2 

' 10.75 
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3/2 



1/4 



(B9) 



9*{ t <j)_ 

By going beyond the sudden decay approximation, one 
finds that AN^ aa is overestimated by about 10%. Ac- 
counting for this by multipling (|B8[) with 0.87, 



AN^ aa 



(T) 



0.82 (\ GeV 



,.2 



1/2 



/PQ 
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Y 



cq/TP ' 



[x 2 + 2(a s /7r) 2 ] 3 / 2 



10" 



(B10) 



(0 


g*s(T)~ 


4/3 g,(T^ 




10.75 _ 


g*s{Ta) 1/3 



This equation together with A = A a ^ 99 A a ^ 99 , as given 
by (|B5p and (|B6p for x > 0.1, can allow for a better 
understanding of the numerical AA^"" results obtained 
in Sect. |IV| 

In Sect. [V] an additional contribution to AN e g in the 
form of axions and axinos from late decays of gravitinos 
is considered. In the sudden decay approximation, 



-™J (|) tfPW/A (BID 



with a temperature at the decay time t = of |20j, |76j, 

iza ' 



which deviates by at most 20% from our numerical results 
for x < 0.02 - shown e.g. by the dot-dashed lines in 



T d = 24 eV 



G 



100 GeV 



3/2 



(B12) 
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Indeed, the gravitino decay happens always after BBN 
so that g*s{T) — g*s(Tg) = 3.91 for the values of mg 
considered in Sect. [V] Proceeding as above, we obtain 



0.42 / 100 GcV 
A V 



1/2 



G 

ht 11 



(B13) 



after multiplying (IB11[) by a factor of 0.87 again to com- 
pensate for overestimation by the sudden decay approx- 
imation. Without this correction and for A = 1, the 
above estimate agrees with the one given in Ref. [20| . 
To understand better the numerical results of Sect. [Vlin 
which AiVeff ~ AN^ aa + AN^ aa , (|BT0)l and (|BT3)) 
can be used together with A = A a ^ 99 , as given by (|B6[) 
for x > 0.1 or by (fB7l) for x < 0.02. 



Appendix C: Description of Changes for m„ 7^ 

Throughout the main part of this work, we use rag — 
m a . In this Appendix we describe the differences with 
respect to the rag — m a case that one faces if those two 
masses differ. Note that we limit our discussion below to 
the regime with r CT < 1 s. In particular, we do not ad- 
dress additional cosmological constraints [18] appearing 
towards small m a that imply longer lifetimes and thereby 
saxion decays during/after BBN. 

Towards very large m a , additional decay channels into 
sparticles may open up such as the decay a — > gg, whose 
width can be derived from (El), 



a 2 m CT m| 



PQ 



2nig 



(CI) 



If this decay occurs after the freeze-out of the LOSP, the 
produced gluinos will not be thermalized. Each of them 
will lead to one LSP and thereby contribute 

n^h 2 = 2ra LSP BR(a -> gg)Y^/ TP (T low )s(T )h 2 /p c , 

(C2) 



similarly as discussed for axino decays in Sect. IIVI In- 
deed, a significant excess over (fT4l) in the G CDM case 
can again be avoided when the saxion decays prior to the 
LOSP decoupling. Because of the additional decay chan- 
nels into axions and gluons, this will be easier to realize 
than for the axino decay in Sect. IIVI In the axion CDM 
case with the very light axino LSP, (|C2j) can be much 



smaller because of a much smaller jtilsp = m a ^ 37 eV. 
Here however such decays could lead to additional con- 
tributions to AiVoff in the form of relativistic axions. If 
the saxion decays prior to the LOSP decoupling, again no 
additional constraints will be expected. A more detailed 
discussion of effects related to the a — > gg decay is left 
for future work. In the following description of changes 
such effects are assumed to be negligible. 

In the G CDM case considered in Sect. IIVI increasing 
(decreasing) m a relative to a fixed value of mg as indi- 
cated on the horizontal axes moves the AN e g contours to 
the left (right) in Figs.[2ja)-(c) and[3l For the x — 1 case 
presented in Figs.[2ta)-(c), there is practically no change 
of the f2g/i 2 contour since the entropy released in saxion 
decays is negligible. For x < 1, A depends on m a as can 
be seen in Fig. [ljb). For a fixed rag value, increasing 
(decreasing) m a reduces (enhances) the dilution due to 
saxion decay and thus results in more (less) restrictive 
upper limits on Tr imposed by f2g/i 2 < 0.129, i.e., the 
respective contours will more downwards (upwards) and 
show a less (more) pronounced dip. 

In the a CDM case considered in Sect. [V] increasing 
(decreasing) m a relative to a fixed rag value as indicated 
on the horizontal axis moves the dotted AiV; 



cir 



con- 



tours in Fig. [|3a) to the left (right). The solid contours 
depicting the sum of both extra radiation components 
change accordingly in Figs. [5ja) andEJa), while there is 

no effect on the dashed AN^[* aa contours for x = 1. 
For x < 1, the entropy release from saxion decays can 
become sizable. For a fixed mg value, increasing (de- 
creasing) m a gives a smaller (larger) A and affects the 
A~N'^ aa contours in a way that is qualitatively compa- 
rable to their change observed for increasing (decreasing) 
x. In turn, the dashed, dotted, and dash-dotted curves in 
Fig. EJb) that present the sum of both AiV e g components 
change accordingly as well as all AA^g contours shown 
in Figs.ObHd). 
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